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o ! , 

i We describe a consistent approach to factorization of scattering amplitudes for exclusive processes 

/"*\ , beyond the leading twist approximation. The method involves the Taylor expansion of the scatter- 

ing amplitude in the momentum space around the dominant light-cone direction and thus naturally 
C*") ' introduces an appropriate set of non-perturbative correlators which encode effects not only of the 

lowest but also of the higher Fock states of the produced particle. The reduction of original set 
of correlators to a set of independent ones is achieved with the help of equations of motion and 
invariance of the scattering amplitude under rotation on the light cone. We compare the proposed 
method with the covariant method formulated in the coordinate space, based on the operator prod- 
uct expansion. We prove the equivalence of two proposed parametrizations of the pr distribution 
O l' amplitudes. As a concrete application, we compute the expressions of the impact factor for the 

transition of virtual photon to transversally polarised p-meson up to the twist 3 accuracy within 
. these two quite different methods and show that they are identical. 

^ I PACS numbers: 12.38.Bx, 13.60.Le 

o : 

I. INTRODUCTION 

O \ 

The study of exclusive reactions in the generalized Bjorken regime has been the scene of significant progresses 
0"\ ' in the recent years, thanks to the factorization properties of the leading twist amplitudes [l[ for deeply virtual 
Compton scattering and deep exclusive meson production. It however turned out that transversally polarized 
p— meson production did not enter the leading twist controllable case Q but only the twist 3 more intricate 
part of the amplitude [!, H, H[. This is due to the fact that the leading twist distribution amplitude (DA) 
of a transversally polarized vector meson is chiral-odd, and hence decouples from hard amplitudes at the 
twist two level, even when another chiral-odd quantity is involved 0] unless in reactions with more than two 
final hadrons [6j. An understanding of the quark-gluon structure of a transversally polarized vector meson 
is however an important task of hadronic physics if one cares about studying confinement dynamics. This 
quark gluon structure may be described by distribution amplitudes which have been discussed in great detail 
0, Q. On the experimental side, a continuous effort has been devoted to the exploration of p-mcson photo 
and electro-production, from moderate to very large energy 0, The kinematical analysis of the final 

7T— meson pair allows then to separate the different helicity amplitudes, hence to measure the transversally 
polarized p— meson production amplitude. Although non-dominant for deep electroproduction, this amplitude 
is by no means negligible at moderately large Q 2 and needs to be understood in terms of QCD. Up to now, 
experimental information comes from electroproduction on a proton or nucleus. Future progress may come 
from real or virtual photon-photon collisions, which may be accessible either at electron-positron colliders or 
in ultraperipheral collisions at hadronic colliders, as recently discussed [TTI. Il2j. 

In the literature there are two approaches to the factorization of the scattering amplitudes in exclusive pro- 
cesses at leading and higher twists. The first approach [1, [l3| , being the generalization of the Ellis-Furmanski- 
Petronzio (EFP) method [l4| to the exclusive processes, deals with the factorization in the momentum space 
around the dominant light-cone direction. We shall call it the Light-Cone Collinear Factorization (LCCF). On 
the other hand, there exists a covariant approach in coordinate space successfully applied in Q for a systematic 
description of distribution amplitudes of hadrons carrying different twists. This approach will be called the 
Covariant Collinear Factorization approach (CCF). Although being quite different and using different distribu- 
tion amplitudes, both approaches can be applied to the description of the same processes. This fact calls for 
verification whether these two descriptions are equivalent and lead to the same physical consequences. This 
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can be clarified by establishing a precise vocabulary between objects appearing in the two approaches and by 
comparing physical results obtained with the help of the two methods. 

The first aim of our paper is to prove that LCCF and CCF are equivalent methods for the description of 
exclusive processes. For that we derive the dictionary between DAs appearing in the LCCF method and in 
the CCF one. We perform our analysis within LCCF method in momentum space and use the invariance of 
the scattering amplitude under rotation of the light-cone vector n tx , which we call n-independence condition. 
This method leads to a definitions of relevant soft correlators which are generally not independent ones. The 
reduction of their number to a minimal set of independent correlators is obtained with the use of equation of 
motions and of the n-independence condition. We obtain the same number of independent correlators in both 
LCCF and in CCF approaches and establish explicit relations between them. 

As a concrete application, the second aim of our paper is to calculate within both methods the impact factor 
7* — * Pt, which is the building block of the description of the 7* p — > pp and 7*7* — > p p processes at large s , 
up to twist 3 accuracy and to verify that we get a full consistency between the two results. 

The paper is organized as follows. In Section |TT] we discuss the general framework of the LCCF method. 
In Subsection III A[ we recall basics on factorization within the LCCF method. In Subsection III Bl we present 
the parametrization of the matrix elements relevant to p-meson production. In Subsection III CI we derive the 
constraint on these matrix elements coming from the QCD equations of motion. In Subsection III Dl we derive 
additional constraints based on the n-independence condition and we then perform the reduction to a minimal 
set of distribution amplitudes. This results in the dictionary given in subsection lll El In Section HTTl we compute 
the 7* — > pt impact factor. After recalling the necessary definitions and kinematics in Subsection IIII Al we 
firstly perform the calculation in the LCCF framework in Subsection IIII Bl , then in the CCF framework in 
Subsection IIII CI We compare the two approaches in subsection IIII Dl Section [IV] presents our conclusions. A 
few appendices present the calculational details needed to complete the proofs. Partial results of this paper 
have been briefly presented in 0, [l6| . 



II. FACTORIZATION OF EXCLUSIVE PROCESSES IN THE LIGHT-CONE COLLINEAR 

APPROACH 



A. Factorization beyond leading twist 



Let us start with the most general form of the exclusive amplitude for the hard process A — > p B (where A 
and B denotes initial and final states in kinematics where a hard scale allows a partonic interpretation) which 
we are interested in, written in the momentum representation and in axial gauge, as 



A= I d 4 ltY 



+ / d% d 4 £ 2 tr 



(1) 



where H and are the coefficient functions with two parton legs and three parton legs, respectively, as 
illustrated in Fig[T] on the example of 7* — -> p impact factor, which is defined in full length in Sec lIIII In 
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Figure 1: 2- and 3-parton correlators attached to a hard scattering amplitude in the example of the 7* — > p impact 
factor, where vertical lines are hard t— channel gluons in the color singlet state. 



(fTj), the soft parts are given by the Fourier-transformed two or three partons correlators which are matrix 
elements of non-local operators. We consider the leading asymptotics of l/Q expansion, separately for the cases 
of longitudinally (twist 2) and transversely polarized (twist 3) meson production. The amplitude ([1]) is not 
factorized yet because the hard and soft parts are related by the four-dimensional integration in the momentum 
space and by the summation over the Dirac indices. 

To factorize the amplitude, we first choose the dominant direction around which we intend to decompose our 
relevant momenta and we Taylor expand the hard part. Let p and n be the conventionally called "plus" and 
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"minus" light-cone vectors, respectively, normalized as p ■ n = 1 . We carry out an expansion of I in the basis 
defined by the p and n light-cone vectors: 

tin = yiPn + (Zi-p)'n>ii+ti~ li , Vi=ti"n, (2) 

and make the following replacement of the integration measure in ((T|) : 

d A l l ^dH t dy l 5{y l -i-n). (3) 

Afterwards, the hard part coefficient function H(£) has to be decomposed around the dominant "plus" direction: 

(t-yp) a + ... (4) 

l=yp 

where (£ — yp) a ~ ta up to twist 3. One can see that the above-mentioned steps ©-(H]) do not yet allow us 
to factorize collinearly the amplitude in the momentum space since the I 1 - dependence of the hard part is an 
excursion out of the collincar framework. To obtain a factorized amplitude, one performs an integration by parts 
to replace by d£ acting on the soft correlator. This leads to new operators O 1 - which contain transverse 
derivatives, such as ^d^tp, and thus to the necessity of considering additional DAs $ (£)■ This procedure 
accomplishes the factorization of the amplitude in momentum space. Factorization in the Dirac space can be 
achieved by the Fierz decomposition. For example, in the case of two fermions, one should project out the Dirac 
matrix tp a (0) ippi 2 ) which appears in the soft part of the amplitude on the relevant T matrices. Thus, after all 




Figure 2: Factorization of 2-parton contributions in the example of the 7* — > p impact factor. 




Figure 3: Factorization of 3-parton contributions in the example of the 7* — > p impact factor. 



these stages, the amplitude takes the simple factorized form 1 , 



A= / dytv[H qq (y)T] <f> r Jy) + / dytv 



H^(y)r\ + J d yi dy 2 ti[H^ g ( yi ,y 2 )T] ^ gfl (yi, y 2 ) , (5) 





in which the two first terms in the r.h.s correspond to the two parton contribution and the last one to the three 
body contribution. This is illustrated symbolically in the example of the 7* — * p impact factor in Fig[2] for 
2-parton contributions and in Fig[3]for 3-parton contributions. 

Alternatively, combining the two last terms together in order to emphasize the fact they both originate from 
the Taylor expansion based on the covariant derivative, this factorization can be written as 



± 1 
A = J dytv[H(y)T] $ r (y) + J dy x dy 2 tr [H»( yi , y 2 ) T] ^(yi,y 2 ). 



(6) 



Despite the fact that these formulae are given here up to twist 3, the method can be extended to higher twist contributions. 
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For definiteness 2 , let us focus on the p- meson production case where the soft parts of the amplitude read 

* r (2/) = J £e-^(p(p)|^(An)r^(0)|0) (7) 

— OG 
+ 00 

*£(l/i,to) = / ^^e-^^-^-y^ipip^Xt^Ti Bj(A 2 n)^(Q)|0). 
—00 

where 

iD^idij.+gA^. (8) 

Eq.([7]) supplemented by the appropriate choice of the Fierz matrices defines the set of non-perturbative corre- 
lators relevant for the description of the p— meson, which we will now discuss. 



B. Parametrizations of vacuum— to— p- meson matrix elements up to twist 3 

In this section, we introduce the parametrizations of the vacuum-to-p-meson matrix elements needed when 
calculating the process of exclusive p~ production. As a concrete example, we shall below calculate the 7^ — > pr 
impact factor. Since we will follow two different approaches for our calculations, it is instructive to present two 
ways for parametrizing the corresponding matrix elements. 



1. LCCF parametrization 

We insist on the fact that in LCCF approach, the coordinates Zi in the parametrizations have to be propor- 
tional to the light-cone vector n. This is in contrast to the CCF approach where z lies on the light cone but 
does not correspond to any fixed light-cone direction. The transverse polarization of the p— meson is defined by 
the conditions (at twist 3, p p ~ p) 

€t ■ n — ct ■ p — , (9) 

i.e. has only a _L component, while has no _L component. Now, we introduce the parametrizations of the 
vacuum-to-yO-meson matrix elements needed for the calculation, for example, of the 7^ — > pt impact factor. 
Keeping all the terms up to the twist-3 order with the axial (light-like) gauge, n ■ A = 0, the matrix elements of 
quark-antiquark non-local operators can be written in terms of the light-cone basis vectors as (here, z = Xn) 

09(^)^(2)7^(0)10) = m p f p [Mv) (e* ■ n) Pfl + <p a ( y ) e* Tp ] , (10) 

(p(Pp) I'0(^)T5Tp'0(O)|O) = m p f p i(p A (y)e fia ( 35 e^ a p l3 n s , (11) 

where the corresponding flavour matrix has been omitted 3 , and where we use e 0123 = — e 012 3 = 1 and 75 = 
i 7 7 1 7 2 7 3 . For the sake of conciseness, we denote = the Fourier transformation with measure 



/ dyexp[iyp- z] , 
Jo 



(12) 



where z = Xn. The momentum fraction y (y = 1 — y) corresponds to the quark (antiquark). Note that the 
decomposition over the 7-matrix basis has been taken in the form: 

- 0>i>) = i(iHV>)7 M + ^757^)7^75 + ■••, (13) 



2 In the following, the notations |p) or | V) will be used when the specific nature of the vector meson does not matter. 

3 The normalization in (1 1 1 1 1 P thus corresponds to a meson which would be a one flavour quark— antiquark state |V) = |//), with 

for example (V / (p)|"i/7(z)7 M i/'/(0)|0) J m v f v [---]. 
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in such a way that the minus sign in front of the axial term is absorbed into the axial correlators. The matrix 
elements of the quark-antiquark operators with transverse derivatives are parametrized according to 

{p(p p )\i>{z) llx i dl 0(0)10) = x m p f p vl(y) Pll e* Ta (14) 
{p{p P )\^{ z )l5l^i dl ip(0)\0) = m p fpiip^iy) e aX p8 e^P 13 n s , (15) 

where we introduced d p = \(d p — d p ) which is the standard antisymmetric derivative. The DAs ipi, 923, pa 
satisfy the normalization conditions 

11 1 
dypi(y) = l,f dyip${y) = l and J dy (y - y) <p A (y) = ^ ■ (16) 


In the same way, the matrix elements of quark-gluon non-local operators can be parametrized as 
(p(PpM(zih P gAl(z 2 )ij(0)\0) = m p f% p B{y l ,y 2 ;y g )p p ,e* Ta , 

(p(p p )m Zl )j 5lp gAl(z 2 )^(0)\0) = m p fi p iD{y l ,y 2 -y g ) Pp e aX ^e* T x pPn\ (17) 

where the momentum fractions y%, y 2 and y g correspond to the quark, antiquark and gluon, respectively. The 
symbol = now stands for (here, Zi — Xi n) 

1 

J dy 1 dy 2 dy g S(y 2 ~ yi - y g )exp[iy 1 p ■ z 1 +iy g p- z 2 ] . (18) 


In the r.h.s. of (|17[) . it is useful to perform the integration over the gluon fraction y g (which then equals y 2 — yi). 
Afterwards, the parametrizations (jTTJ) take the forms: 

( P (p p m Zl )j p gAl(z 2 )^(0)\0) = m p f? p B( yi ,y 2 ) Ptl e* Ta , (19) 
(p(p p )\i>( Zl ) l5lp gAl(z 2 )i;(0)\0) ^ m p tf p W(y u y 2 )p p s aXf3S e* T x p^n 5 , (20) 

with the symbol = implying 

1 1 

dyi J dy 2 exp[iyip- z 1 +i(y 2 -y 1 )p- z 2 ] . (21) 


Note that the positivity of the gluon light-cone momentum fraction imposes that quark-gluon parameterizing 
functions have the form 

def def 

B(yi,y 2 ) = B(yi,y 2 ;y 2 - yi)9{y x < y 2 < 1), D(y 1: y 2 ) = V(y 1 , y 2 ; y 2 - yi) 9(y 1 < y 2 < 1) . (22) 
As we already mentioned by writing Eq.([5]), it is also natural to introduce the following objects: 

(p(p P M( z ih^ D a (za)i>(0)\0) = m p } P B{y u y 2 )p p e* Ta , 

(p(p P M (zi)l5l p,i Dl (z 2 )-0(O)|O) = m p f p it)(y ll y 2 )p p e a xi38eT*p fs n s 1 (23) 
where these parameterizing functions are now equal to 

B{yi,y 2 ) = ^(^pj {yi) + <fj (y2)^j 6(y! - y 2 ) + B( yi ,y 2 ), 

j/a) = \ ^a(j/i) + ¥>5(w)) %i - Vi) + (3 D(y u y 2 ) , (24) 
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with the dimcnsionless coupling constants 

C v = !Li and C 3 A = ^- (25) 

Note that the function ipi corresponds to the twist-2, and functions B and D to the genuine (dynamical) twist-3, 
while functions tps, tpA,tpf, tp\ (or alternatively B and D) contain both parts: kinematical (a la Wandzura- 
Wilczek, noted WW) twist-3 and genuine (dynamical) twist-3. 

In (|10 [) — (|li p . the functions (p-y, <p$, ifA, f'l and Lp\ parameterizing the two-particle correlators obey the 
following symmetry properties: 

<pi{y) = <pi(l- y), <P3(v) = ^3(1 - y), va{v) = -<pa(i - y), <pT(y) = - v), <Pa(v) = - v) ■ (26) 

These symmetry properties result from G-conjugation (or C-conjugation for neutral mesons). At the same time, 
the symmetry properties of the functions parameterizing the quark-gluon correlators are: 

B{yi,V2\Vg) = -B(l-J/2, 1 -VDVg), v(vi,V2;Vg)='DQ--V2,i-vi;vg)- (27) 

Notice that, in the case of three-particle functions, G-conjugation involves the replacement: y\ <-*• j/2, while the 
gluon fraction y g remains invariant under G-conjugation. 



2. CCF parametrization 

We recall and rewrite (doing standard fields transformations) the original CCF parametrizations of the p 
DAs 7], adapting them to our case when vector meson is produced in the final state. The formula for the 
axial- vector correlator reads 



(p(Pp)H( z ) [z, 017^75^(0)10) = ^fprnpEf, 



e T p z 



(28) 



where we denote 4 



e T pz afij * 



(29) 



and where 



[z\, z 2 ] = Pexp 



ig J dt (zi - z 2 )„A>*(t zi + (1 - t) z 2 ) 




(30) 



is the Wilson line, defined in accordance with the convention ([8]). The transverse vector ex is orthogonal to the 
light-cone vectors p and z. Neglecting mass effects, i.e. up to twist 3 level, it is decomposed as follows 



erfj, 



e ■ z 



p-z 



e • p 
p ■ z 



(31) 



where e is the meson polarization vector. Since as we will discuss later n can be arbitrary, and since the concrete 
definition of n influences the definition of transverse polarization, it is useful to remove the dependence on ct 
in correlation functions. For that, we use Eq. (|3ip and rewrite the original CCF parametrization in terms of 
the full meson polarization vector e. This is already done for the axial-vector correlator (|28| since due to the 
properties of fully antisymmetric tensor e^ap one can use in the r.h.s. of (|28[) the full meson polarization vector 
e instead of eT- 



4 Note that, as already emphasized, our sign convention for the antisymmetric tensor is e 0123 = 1, opposite to the one used in 
Ref. 0. The corresponding sign change is taken here into account. 
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The definition of 2-parton vector correlator of a p-meson can be written in the form 

1 



(p(jp P )\ip(z) [z, 0] 7^(0)|0) =f p m p J dye 





iy{p-z) 



P„ 



h (y) + e* Tt . g { ± (y) j^r^ 



9 * 

m e ■ z 



p ■ z 



93 (y) 



(32) 



All distribution amplitudes describing two particle correlators are normalized to unity 



dy [4>\\,g±\g±\g3^ (y) 



= l 



(33) 



Using relations (|311 133|) and integration by parts one can rewrite the vector correlator (132p in the form 



(p(p P M(z) [z, 0] 7^(0)|0) = f p m p / dye*«<r* 



-ip p (e* ■z)h(y) + e* fJ/ g i l'{y)+iz f _ l 



2 p- z 



My) 



(34) 



where we introduce the auxiliary functions 



h(y) = J dv (4>\\(v)-gW(v) 
o 

y 

h(y) = J dv {gs(v)-g^\v)) 



(35) 



(36) 



Note that the r.h.s of (|34[) now only involves the full polarization vector e, as was noted above for the axial 
correlator. The last term in the r.h.s. of (|3"4"|) contributes to the physical amplitude starting from the twist 4 
level only, therefore we will neglect it in the following. 

For quark-antiquark-gluon correlators (up to twist 3 level) the parametrizations of Ref. [7j have the forms 5 



(p(p P )\i>( z )[z,t z]j a g G pv (t z)[t z,0]V>(0)|0) = -ip a \p M e^ u - p„e* ±IJ ]mpf 3p / DaV(a 1 ,a 2 )e l(p ' 



i(p-z)(a 1 +t a.g) 



, (37) 



(p(p p )|^(z)[z,tz]7a75ffG^(iz)[^,0]^(0)|0) = -p Q b P el,-p, e lJm,/ 3 A „ / Da A(a u a 2 )e l ^ z ^ +t a ^ , (38) 



P \" l p iz P 

where a±, a 2 , a g correspond to momentum fractions of quark, antiquark and gluon respectively inside the 
p— meson, 

ill 

C-S0) 

ooo 



iii 

J Da — J da\ J da 2 J da g 8(1 — a\ — a 2 — a g 



and G M „ = — \e pV upG a P . These three partons DAs are normalized as follows 



J D a(a 1 -a 2 )V(a 1 ,a 2 ) = l, J Da A(a x , a 2 ) = 1 



(40) 



In what follows we will work in the axial gauge A ■ n = 0, n 2 = 0. In this gauge the gluon field can be expressed 
in terms of field strength as follows 



A a (y) = J dae" e,7 n' 3 G a 0(y + an) 



(41) 



5 Note that in those definition f p , fY and f£- have dimension of mass. This is agreement with Ref. Q but differ from Ref. 0] in 
which fY p and f^ p have dimension of mass square. 
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which implies that the (q A q ) correlators involving the gluon field A reads 

(p(p p )\iP(z)^gA a (tz)m\Q)=-P»e* Ta m p f? p f — e *Cf»X«+*«.) V{cc x , a a ) , 



2pe T 



Da 



(42) 
(43) 



C. Equations of motion 

The correlators introduced above are not independent, since they are constrained by the QCD equations of 
motion for the field operators entering them (see, for example, j4;])- In the simplest case of fermionic fields, they 
follow from the vanishing of matrix elements ((iD(0)i/j(0)) a 'ipp(z)) = and (ip a (0)i(D(z)ip(z))p) = due to 
the Dirac equation, then projected on different Fierz structure. 

Let us start with the QCD equation of motion written for the fermion field ip(0): 



(ip (0)V(0) 0, 



(44) 



where (. . .) denote arbitrary hadron states which we here specify as (p\ . . . |0). Also, we stress that, in (|4"4"|) . 
the fermion fields ip and tp should be understood as fields with free Dirac indices. Then, we first focus on the 
quark-antiquark part of (]44[) which can be written as 



(45) 



Here, we separate out the longitudinal derivatives from the transverse ones. Working with the longitudinal 
derivative contribution, we get 



. j4 —iyp-z — iyp-x 

i"f p laze ' 



QL 

dx, 



■(tI>{x)$[z)) 



y$ / d ze 



iyp-z — typ-x 



(46) 



where an integration by parts has been used. Let us now decompose (ip(x) ip{z)) over the 7-basis (the Fierz 
decomposition) : 



(ip(x)ip(z)) = -{i()(z)j a ip(x))j a + ^(ip(z)-/5l a tp(x))-/ a l5 ■ 



(47) 



With (|4"7) . after the use of the parametrization of the relevant correlators, one gets for the longitudinal derivative 
contribution (see, (|46|) ): 



|iM T v?3(y) + i-i>f ^I5Va{v) = --f pe ?\ y<pa{y) + y<PA(y) 



(48) 



where again we introduced the short-hand notations: 



®p e T — ®ol (3 V e T ' ^Tp — Zpe^pn 

We thus have for the longitudinal derivative contribution: 



(49) 



d i ze -iVP-z-iyp-x q q ^x)^(z)) 



x=0 



--v peT m p f p \y <p 3 (y) + y<p A (y) 



(50) 



While, the correlators with the transverse derivatives in (|45p can directly be expressed via the corresponding 
parameterizing functions with the help of (Ti"4"|) and (fT5]) : 



-iyp-z-iyp- 



*{l $ T TP(x)ij(z)) 



x=0 



= -^ P e* T m p f p <{ ip x (y) + (f A (y) 



(51) 



Therefore, within the WW approximation (where all genuine twist 3 are disappeared) the equations of motion 
takes the following simple form: 

vl(y) = -y<p¥ w (y), (52) 
where the plus (minus)-combination is defined as 

¥>±(y) = ^"vector" (y) ± P "axial" (v)- (53) 

Let us now take into account the quark-gluon correlators. Using (|17p . one can obtain that 

- lp (A T P(0)m^)) = ll P i^(z)laA T P(0)i:(0))j a + ^(z) l5la A T P(0)mhaJ 5 }, (54) 



where 

1 

V = " l pJ3p Pa^Tp 



^{z) la gA T p {Q)i>{Q))=m p g pPa e* Tp / dy x dy 2 e^ z B( Vl , y 2 ), 



o 

l 



l 

dx 



(i>{z) l5la gA p {0)i>{0))=m p f£ p i Pa a Tp J d yi dy 2 e^P z D( yi ,y 2 ). (55) 

o 

Thus, combining the quark-antiquark and quark-gluon correlators, one derives the following relation: 

(B{y,x)+D(y,xj)=-yip + {y), (56) 
o 

where one used the notations ([24)1 . 

In a similar way, we can derive the relation associated with the C-conjugated equation: 

(i/)(0)$(z)if) (0)> = 0. (57) 

which reads 

l 



dx 



(B(x,y)-D(x,yj)=ytp-(3,). (58) 







Combining the relations (|5"6"|) and (|58p with the use of (f2"4"|) and (JHSJ), we obtain 

m *P3(yi) + m pa(vi) + <pT(yi) + fKvi) 



i 



= -J d yi i$ B (y^ + # £>(w, y 2 )] (59) 



and 

yi ^3(2/1) - yi VA{yx) - <pf{yi) + <pa(vi) 

dy 2 [-C 3 y B(y 2 , Vl ) + ^ D(y 2 , Vl )] . (60) 



Note that Eq. (fBT)]) can be obtained by the replacement y\ — > y\ in and the use of symmetry properties (|26[ 
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D. Additional set of equations 

1. Light-cone factorization direction arbitrariness 

Contrarily to the light-cone vector p related to the out-going meson momentum, the second light-cone vector 
n (with p ■ n = 1), required for the parametrization of the needed correlators introduced in section Hi B II is 
arbitrary. The physical observables do not depend on the specific choice of n, thus the scattering amplitudes 
should be n— independent. For any specific process, there is a natural choice for n, which one may denote as 
riQ. For instance, in forward e — p collision, the proton momentum defines no- More generally, one may expand 
an arbitrary choice of n as [U, [lj] 



with the two constraints 



p ■ n 



tp f _ l + (3n 0fl , + n^ , (61) 



1 and n 2 = , (62) 



which fixes the coefficients /3 = 1 and a = — . The light-cone vector n is thus parametrized by its transverse 
components n±_ . 

Let us now analyse the various source of n— dependence. First, it enters the definition of the non-local 
correlators introduced in Sec lIIB II through the light-like separation z = Xn. These correlators are defined in 
the axial light-like gauge n ■ A = , which allows to get rid of Wilson lines. Second, it determines the notion 
of transverse polarization of the p . Last, n inters the Sudakov decomposition |(5J) which defines the transverse 
parton momentum involved in the collinear factorization. Note that this notion of parton transverse momentum 
should not be confused with the notion of transverse momenta of external particles (e.g. in the case of 7* — > p 
impact factor to be discussed in Sec lIIIl entering in 7* N(p 2 ) — > p(p) N process, the t— channel gluons have a 
transverse momentum determined within another Sudakov basis defined by the external light-cone momenta 
Pi = p and incoming nucleon momentum p 2 ). 

This n— independence principle leads to additional non-trivial constraints between the non-perturbative cor- 
relators entering the factorized amplitude. It was crucial for the understanding of inclusive structure functions 
properties at the twist three level [14j and its relevance for some exclusive processes was pointed out in [4| . We 
show now that this condition expressed at the level of the full amplitude of any process can be reduced to a set 
of conditions involving only the soft correlators. The obtained equations are process independent and do not 
assume a priori any Wandzura-Wilczek approximation. The strategy for deriving these equations relies on the 
power of the Ward identities to relate firstly amplitudes with different number of legs and secondly higher order 
coefficients in the Taylor expansion (T5]) to lower order ones. 

In the case of processes involving px production up to twist 3 level, we will now derive the equations 

1 

d . Tf s , s ( , tV f dy 2 



dyi J yi-y\ 



d f 

^"^(yi) - <pa(vi) + J 





1 

dy 2 



V2 - yi 



(D(y u y 2 )+D(y 2 , yi ))=0. (64) 



The n— independence of A for an arbitrary fixed polarization vector e is expressed by the condition 

d 



dn^ 



A = 0, (65) 



which we write in the form 



d . dn a dA d(e*-n) dA . a a . d . , 8A 

d^ A = ^^ + -^dW^) = [ ~ n ^ p +9 ^ ] d^ A+e ^dW^)' (66) 

Let us emphasize the fact that although n fixes the gauge, the hard part does not depend on this gauge fixing 
vector, as we will show below after the technical derivation of Eqs. ([63ll64|) . It means that the variation of n only 
affects the n— dependence related to the definition of transverse momentum £± and of transverse p polarization. 
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We also note that the appearance of the total derivative in Eqs. (j6"5l may be interpreted as a (vector) 
analog of the renormalization group (RG) invariance equation when the dependency on the renormalization 
parameter coming from various sources cancel. One can view this as a RG-like flow in the space of light-cone 
directions of contributions to the amplitude where the polarization vector plays the role of a beta function. 

The scattering amplitude A receives contributions from the vector correlators, which result into j[ vector ; and 
from the axial vector correlators, which lead to J[ axtal part of A. Due to different parity properties of the vector 
and the axial-vector correlators, the condition ([55]) means effectively two separate conditions: 



d 



dn± 

and 

d 



dn. 



A vector = q ( 6? ) 



The dependence of A on the vector n± is obtained through the dependence of A on the full vector n. This 
dependence on n is different in j[ axtal and in j[ vector parts. 

The dependence of j[ axtal on the vector n enters only through the expression e pn ^ 1 involving the contraction 
with the momentum p and in which the indices f3 and 7 are contracted with some other vectors. Thus the 
condition ([55]) is equivalent to 

dn^_ 8 A<xxlal = + A axial = ^rA axial = (69) 

on^_ on a on a 0717 

where we took into account the peculiar dependence of j[ axial on n discussed above. This will lead at the level 
of DAs to the equation ([64]l . 

In the vector part, the dependence with respect to n is identified by rewriting the polarization vector for 
transversally polarized p with the help of the identity 

e^, = e^-p^e-n, (70) 

since e • p = . Thus the dependence of A vector on the vector n enters only through the scalar product e* • n, 
and Eq. (j67|) can be written as 

Avector * fi Avector r\ /rri \ 

which results in 

■ A vector = (72) 



d{e* 



from which follows Eq. (p5)l . 

We will now derive equations (|63l [64f using as a tool the explicit example of the 7* — > p impact factor. We 
want here to insist on the fact that the proof is independent of the specific process under consideration, and only 
rely on general arguments based on Ward identities. For the 7* — > p impact factor, which will be computed in 
details in section HTIl one needs to consider 2-parton contributions both without transverse derivative (illustrated 
by diagrams of Figf4]) and with transverse derivative (see FigfS]), as well as 3-parton contributions (see FigslH 
[3 [H]) . Note that these drawing implicitly assume that the two right-hand side spinor lines are closed on the the 
two possible Fierz structures j) or fy* involved in the correlators of p— meson DAs. 

In the color space, each of those diagrams can be projected in two parts, characterized by the two Casimir in- 
variants Cf and N c . The equations (|63l I64|) are obtained by considering the consequence of the n— independence 
on the contribution to the Cp color structure. The n— independence condition applied to the iV c structure is 
automatically satisfied and does not lead to new constraints, as we show in Appendix [AJ 

We start with the derivation of Eq. (|63j) . which corresponds to the vector correlator contributions with Cf 
invariant. The 3-parton (qqg) contribution and the 2-parton contribution involving $ x to A can be reduced to 
the convolution of the leading order hard 2-parton contributions with linear combination of correlators, thanks 
to the use of the Ward identity. 
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(e) (f) 

Figure 4: The 6 hard diagrams attached to the 2-parton correlators, which contribute to the 7* — > p impact factor, with 
momentum flux of external line, along pi direction. These drawing implicitly assume that the two right-hand side spinor 
lines are closed on the the two possible Fierz structures j) or j> 7 s . 

In the case of the 3-parton vector correlator (|19[) . due to (|7D|) the dependency on n enters linearly and only 
through the scalar product e* - n . Thus, the action on the amplitude of the derivative d/dn± involved in (|T2[) can 
be extracted by the replacement e* — ► — p a , which means in practice that the Feynman rule (using conventions 
of Ref. [l7| for computing the T matrix element) gt a "f a e* entering the coupling of the gluon inside the hard 
part should be replaced by — gt a 7° p a . Then, using the Ward identity for the hard part, it reads 

(yi - 2/2)tr [H% 9g (y 1 ,y 2 )p p $] = tr [H^yx) f[ - tr [H q g(y 2 )p] , 

which can be seen graphically as 




as we will show below with more details, and will give the last term of r.h.s of Eq.Q63|). The proof can be settled 
easily relying on a graphical rule in order to use the collinear Ward identity. Indeed within the conventions of 
fl7j |. the collinear Ward identity can be symbolically written as 



1 



{V2-Vi)p (V2-Vi)p 
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where each fermionic line is a propagator. The wavy lines with double arrows are there in order to fulfill 
momentum conservation, since the incoming momentum is y\ p while the outgoing momentum is p ■ 

Let us consider the 3-parton "Abelian" diagrams, illustrated in FigEl Applying the Ward identity (|74|) to 



14 





Figure 6: The 12 "Abelian" (i.e. without triple gluon vertex) type contributions from the hard scattering amplitude 
attached to the 3-parton correlators for the 7* — > p impact factor, with momentum flux of external line, along pi 
direction. 
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(etGl) (etG2) (ftGl) (ftG2) 

Figure 7: The 12 "non-Abelian" - (with one triple gluon vertex) contributions from the hard scattering amplitude attached 
to the 3-parton correlators, for the 7* — ► p impact factor, with momentum flux of external line, along pi direction. 




(gttGl) (gttG2) (httGl) (httG2) 



Figure 8: The 4 "non-Abelian" -(with two triple gluon vertices) contributions from the hard scattering amplitude attached 
to the 3-parton correlators, for the 7* — > p impact factor, with momentum flux of external line, along pi direction. 
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graphs (aGl) and (aG2) gives 




(75) 



where the indicated momentum fractions correspond to flow along the momentum pi of the p~ meson. The last 
line of Eq. ([75)) has been obtained after cancellation of the first and fourth term in the second equality, and the 
two remaining diagrams have been relabelled (in the first term of the last line, one does y% — 1 — ► y% — 1 for the 
outgoing antiquark, and in the second term, one does y\ — > yi for the outgoing quark), as far as the external 
lines are concerned, after using the fact that this does not change its internal structure. 

The same identity applies for each couple of graphs (bGl, bG2), (cGl, cG2), (dGl, dG2), (eGl, eG2) and 
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(fGl, fG2). This leads to the following identity 



dyi J dy 2 B(yi,y 2 ) 



o o 



Pp. 



V2 - VI 

V2 — 1 



+ ^ 




dyi / dy 2 



o o 



B{y\,V2) 
V2 - yi 



l/l—l 



- (yi <-> 2/2) 



= / dy x 



dy 2 



2/2 - 2/i 



[£(3/1,3/2) + S(y 2 ,2/i)] 



+ 



yi - 1 



TT 



TT 



(76) 



where the last line is obtained after performing the change of variable y\ <-» y 2 in the second term. Note that 
the hard part given by diagrams inside Eq. (f?6")) is convoluted with the last term of Eq. (|6li)) . 

A similar treatment of 2-parton correlators with transverse derivative whose contributions can be viewed as 
3-parton processes with vanishing gluon momentum leads to 



dyi / dy 2 S(y 1 - y 2 ) V>1 {Vl) 



Pp. 



V2 - VI 

V2 — 1 



iCTT 



= / dy x / dy 2 8{yi - y 2 ) 



V2 - yi 
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1 




-(77) 



where the last line is obtained after integration by part. This leads to the convolution of the first term of the 
l.h.s of Eq.([63]) with the [• • • ] part in (J77|). 

The second term, with ip\ , of the l.h.s of Eq. (|63[) originates from the 2-parton vector correlator and corresponds 
to the contribution for the longitudinally polarized p with ~ p. The third term with if 3 corresponds to the 
contribution of the same correlator for the polarization vector of pt written as in Eg. ([70)) . To get Eq. ([55|) . we 
used the fact that each individual term obtained above when expressing the n— independence condition involve 
the same 2-parton hard part, convoluted with the Eq. (|63[) through an integration over y\ . The arguments used 
above, based on the collinear Ward identity, are clearly independent of the detailed structure of this resulting 
2-parton hard part. Therefore, we deduce from this that Eq. (J53J) itself should be satisfied. 

A similar treatment for axial correlators leads to Eq.((64|). To prove this, we start from Eq. ([69|) and we note 
that the parametrizations of matrix elements of correlators with axial-vector currents (| 1 51 1171 120p involve the 
quantity 



fj> _a e T p n 



(78) 



in which the index a is contracted with the matrix 7" appearing in the vertex of gluon emission in the hard part 
and the momentum p^ is contracted with the Fierz matrix 7^7 5 corresponding in the hard part to the meson 
vertex. First let us note that in the expression (|78|) one can replace by the full polarization vector e*, i.e. 



I± _a e T p 71 



fi _a e p n 



(79) 



Secondly, the inspection of the quantity (|75|) or ([79")) leads to the conclusion that in order to use the Ward 
identities in a similar way as it was done in the vector part we need to interchange in (|78p the indices fj, *-y a. 
It is done with the help of the Schouten identity, which for our peculiar case means that 



_Q e T pn 



p^e 



fi e T p n 



(80) 



After that, the momentum p a acts on the gluon vertex in the hard part, so the consequences of the 
n— independence of the axial part of the impact factor maybe derived in exactly the same way as we did above 
in the case of vector correlators, since the vector £^e T P n is completely factorized. One then obtains Eq. ([M|) 
from Eq. (|63p after the replacements B — ► D, 923 — > cpA, Lp\ — ► tp T A and (fi — ► since there is no counterpart of 
the twist 2 DA ipi for the axial part. 



Since we rely on the n— independence of the amplitude, one may wonder about the effect of the gauge choice, 
which is fixed by n , on the hard part. The QCD Ward identities require the vanishing of the amplitude in which 
polarization vector of a gluon is replaced by its momentum provided all other partons are on the mass shell. In 
the framework of the kr— factorization (see Sec iIII|) . the t— channel gluons are off the mass-shell. Therefore the 
replacement of the s— channel gluon polarization vector by its momentum leads to the vanishing of scattering 
amplitude up to terms proportional to k\j 1 s where k± are transverse momenta of t— channel gluons. From the 
point of view of the t— channel, the gauge invariance of the impact-factor means that it should vanish when the 
transverse momentum of any t— channel gluon vanishes. To achieve this property it is necessary to include in a 
consistent way not only DAs with lowest Fock state containing only quarks but also those involving quarks and 
gluon, as we will show in detail in Sec lIIIl 

In practice, we here check this invariance by contracting the s— channel emitted gluon vertex in the hard part 
with the momentum, which in collinear factorization is proportional to the p— meson momentum, which leads 
to simplifications in the use of (collinear) Ward identities. 

In order to prove this, one should first project on the various color Casimir structure. In the case of the 
impact factor (see Sec JIIIj) . this means to distinguish N c and Cf terms. In this case, Cf terms arise from 
2-partons diagrams and from 3-partons diagrams where the emitted gluon is attached to a quark line, while N c 
terms are obtained from 3-partons diagrams where the emitted gluon is attached to a quark line only between 
2 t— channel exchanged gluons or from diagrams involving at least one triple gluon vertex. 
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The method is almost identical with the one used above when deriving the n— independence equations. 
Consider first the case of hard 3-partons diagrams entering the vector part of Fierz decomposition. We contract 
the s— channel emitted gluon vertex with its momentum, which is proportional to . The next step is to use 
the same method as the one used in Eq. l|75|) . except that the DA B is not involved here. One thus finally gets 
two groups of 6 diagrams (which differ by the labeling of outgoing quarks, one being y\ and 2/1 for the quark 
and antiquark, respectively, and the other one being 2/2 and 2/2 for the quark and antiquark respectively) . Since 
we started here from the consideration of the — > px transition, each of these 6 hard contributions, due to 
the appearance of the remaining j> from the Fierz structure, thus now encodes the hard part of the transition 
7t ~~ * Pl ■ This transition vanishes in our kinematics, which leads to the conclusion that this hard part is gauge 
invariant. The same treatment can be applied to the hard diagrams with derivative insertion displayed in FigfSJ 
since this insertion corresponds to the peculiar limit of vanishing "gluon" momentum. 

The proof for the axial part of the Fierz decomposition goes along the same line. The only difference lays 
on the appearance of the ]/>j 5 structure, which corresponds to a meson b\ with quantum numbers J PC = l" 1 
instead of 1 for p, leading finally to the hard part of the transition j T — > bi l which vanishes in our kinematics. 

In the case of contribution proportional to N c , one can prove that these hard terms are also gauge invariant. 
This is proven in Appendix [A] The reason is the same as the one which led to the conclusion that N c terms do 
not lead to additional n— independence condition. 

Although our implementation of factorization and n— independence condition is illustrated here on the par- 
ticular example of the impact factor at twist 3, we expect that this procedure is more general and that the 
above method can be applied for other exclusive processes, for which the key tool is still the collinear Ward 
identity. This means in particular that each building block (soft and hard part, for each structure which lead 
to the introduction of a DA) are separately gauge invariant. This fact simplifies dramatically the use of the 
n— independence principle. 

2. A minimal set of non-perturbative correlators 

We now solve the previous equations, namely the two equation of motions (|59I60|) and the two equations 
(|63I64[) coming from the n— independence. This effectively reduces the set of 7 DAs to the set of 3 independent 
DAs <p-L, B, D. 

To start with we represent the distributions <p${y), <Pa(v)i <pT(.V) an d < Va{v) generically denoted as <p(y) as 
the sums 

<p(v) = y ww {y) + v Ben {y) , <p(v) = Mv), Mv), vliv), fl(y) > (si) 

where ip ww (y) and ip gen (y) are contributions in the so called Wandzura-Wilczek approximation and the genuine 
twist-3 contributions, respectively. 

The Wandzura-Wilczek contributions are solutions of Eqs. (|59l [60l [63l [64]) with vanishing 3-parton distribu- 
tions i?(yi,?/2) and D(2/i,2/2), i.e. which satisfy the equations 

m v? w {yi) + m vJT (w) + ww (vi) + <& ww (vi) - o (82) 
yi <pT w (vi) - vi (2/1) - ^ ww (yi) + ww (vi) = ■ (83) 

ww (yi) - -Mvi) + vT^dft) , ww ^ = va w (vi) ■ m 

dyi dyi 

By adding and subtracting Eqs. (|82l83p together with the use of Eqs. ([54]) one obtains equations which involve 
only ipf w and ^ w 

^-<pY w ( yi ) = -(yi - yx)^J w {yx) , 2^) = ^J w (yi) + (m - yi)^<pY w (yi) (85) 

dyi dyi dyi dyi 

and which solutions, satisfying the normalization conditions 

1 1 
Jdy^ w {y) = l and jdyyj w {y) = 1 , (86) 
o o 
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read 



Vl 1 

dv . . f dv . . 

—<Pi{v) ~ / —<Pi(v) 

V I V 



111 



^ w {yi) 



These expressions and Eqs 
1 



Vl i 

dv . . f dv . . 

v J v 

!.n 

T WW A ,J WW 



(87) 



give finally the remaining solutions tp A and tp{ 



, T WW ( \ 

<Pa (2/1 J 



</i / ^yi(w) - yi [ —(fi(v) 

V J v 

yi 



vl ww (yi) 



dv . , f dv 

-yi I —ipi{v) + yx / —ipi{v) 
v J v 

Vl 



We note that these two WW results ((87)) were obtained in Ref. [18| when considering the transition form factor 
B — > p 7 . The distributions ip 9en carrying the genuine twist-3 contributions satisfy the equations 



yi<P 9 3 en (yi) + yi<pT(yi) + rf 9en (yi) + <PA 9en (yi) = - I dy 2 [c^ B( yi ,y 2 ) + ^D( yi ,y 2 )} , (89) 



T gen 



yi ^ 9 3 en (yi) ~ yi ¥>T (Vi) - £ gen (yi) + Va 9en (yi) = - dy 2 B (y 2) Vl ) + C3 4 D(y 2 , y x )] , (90) 



1 

T qen / \ qen / \ /-V f 

^-Vi m) = m) - C 3 J 



1 

d T qen / \ qen / \ /-A f 

-fa^A (yi) = <p A J 



dy 2 



yi - yi 



(B(y 1 ,y 2 ) + B(y 2 ,yi)) , 



dy 2 



2/2 - 2/1 



(D( Vl ,y 2 ) +D(y 2 , yi )) 



(91) 



Similarly as in the WW case, one can obtain equations without distributions tp T gen by adding and subtracting 
Eqs. (Ell E3) together with the use of §2$ 



^ qen / \ , / — , d q en . . 

-i—n (2/1) + {yi-yi)-i— ( PA 2/1 

dyi dyi 
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4C3 4 / -^D^( yil y 2 )-2Cr-^- I dy 2 B(-\ yi ,y 2 )-2Ci^- f dy 2 D^( yi ,y 2 ) , (92) 
J 2/2 - 2/i cm/i J «2/i J 



^ gen / \ . / _ . Q fin , . 

-7— <T A {yi) + {yi - 2/iJ-T— n (2/1 
dyi dyi 



X 11 

H *J y^ Bi+){yi ' y2) ~ 2C *Wi I ^ B(+) (^'^)-2C 3 4 ^ / dy 2 Di-\ yi ,y 2 ) , 



(93) 



where (yi,y 2 ) — 0(yi,y 2 ) ± 0(y 2 ,yi) for O = B, D. Let us note that Eq. (|93|) can be obtained from 
Eq. (|i?2")) and by the interchange <p^ e " <-> (fiff 1 and B ^ D. 

From Eqs. (|92l [93)) supplemented by the boundary conditions B(y, y) = = £>(?/, y) (see later in Section Hi E[) 
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one gets in a straightforward although somehow tedious way the equation for ipf. 



gen 



dy 



qen / \ 



1(1 + 1 

2 \yi yi 



l yi / l yi 

yi J dy 2 -^B(y 1 ,y 2 ) - yi J -^-B(y 2 ,y 1 ) + (yi - yi) I J dy 2 ~~~ + J dv '> 



£(2/2,2/1) 
2/2 - yi 



Hi 



yi / dy 2 —D(yi,y 2 )+yi I —D(y 2 ,y 1 )- / dy 2 



dy 



dyi 



2/2 - 2/i 



D{y2,yi) 

2/2 - 2/i 



.yi 



(94) 



which properly normalized solution can be written in the form 



n kv) = o 



dyi , M2/1 I r 1 
2/1 J m ' 



(95) 



in which both integrals act on the expression inside {...} on the r.h.s of Eq. (|M]) . The expression (|95")1 can be 
simplified after changing the order of the nested integrals. After performing this task we obtain that 



gen / \ 

n (2/) 



(96) 



1 U 1 U 

~f-\f dy^B - &D)(y 2 , u) - f -^(CYB - tfD)(u, y 2 ) - f - - &D)(y 2 , u) 

2 J u [J du J y 2 - u J y 2 - u 



y 
yi 1 



■5 f d 4 I dy 2 ±((ZB + tfD)(u >y2 )- f B + &D)(u, y 2 ) 

2 J u [J du J y 2 - u 



dy 



2 f/ y 



u 



2/2 -u 



Finally, the solution for (pj gen is obtained from the first Eq. (|91[) and 



iPi 9en (y)= I duipf n {u) - (3 I dyx I dy- 



B(yi, 2/2) 

» 

2/2-2/1 



(97) 



y 

The corresponding expressions for (p 9 ^ n (y) and (p^ 9en (y) are obtained from Eq. (f9"6"]) and (f9"7| by the substitutions: 



gen / \ C3 ^ C3 ^ qen / \ 
^3 (2/) < * (2/): 

Tgen/ \ CT^^C^- T 



J gen / \ 
¥>1 (2/) 



_z gen / \ 

<pa (y)- 



(98) 
(99) 



In conclusion of this section, we explicitly succeeded in representing our results l|57|) , l|88p , (IM)) , (|97p , 11980 (|9"9"| 
in terms of 3 independent DAs: the twist 2 DA (p% and the twist 3 DAs B ,D . 



E. Dictionary 

For comparison of expressions (TT91 I2TJ)) with the definitions (T4"2l |4"3"|) we perform the change of variables — ► Z\ , 
— > Z2, ad —> 2/i an d a u = 1 — y 2 , i.e. a g = y 2 — y%. It results in the following identification of the 3-parton 
DAs in LCCF and CCF approaches 

B(yi, y 2 ) = 100 

2/2 - 2/i 

rv \ -4(2/1,1-2/2) , 1fvn 

D{yi,y 2 ) = . (101) 

2/2 - 2/1 
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From (|10Q[ I101[) and Ref.Q follows the boundary conditions B(y, y) = = D(y, y) . 

Taking in Eqs. (j32l [34|) the coordinate z along the light-cone vector n, z — An, permits the identification of 
the vector DAs in Eq. lTUl) : 

Mv) =h(y)> fs(y) = 9 ( lHy), (102) 



and of the axial DA in Eq. (fTT|) 



We checked also the validity of Eqs. p02[ 1 1 03[) directly by the use of our explicit solutions ([571 [Ml and 

expressions for g^ and g^ given Ref. Q in terms of (fiu , V and A DAs. This non trivial check can be done 
with the help of methods similar to those used in Appendix [B] when comparing the results of calculation of the 
7* — > Pt impact factor in LCCF and CCF approaches. 

III. 7* -> P t IMPACT FACTOR UP TO TWIST THREE ACCURACY 

A. General recall on impact factor representation and kinematics 

The 7* — > p impact factor enters the description of high energy reactions in the kx factorization approach. 
As an example, one may study the reactions 

i*(q)+i*(q')^PT(Pi) + p(P2) (104) 

or 

-t*(q)+N ^ p T ( Pl ) + N (105) 

where the virtual photons carry large squared momenta q 2 = —Q 2 (q 12 — ~Q' 2 ) ^qcd > an d the Mandelstam 
variable s obeys the condition s 3> Q 2 ,Q' 2 ,—t ~ r 2 . The hard scale which justifies the applicability of 
perturbative QCD is set by Q 2 and Q' 2 and/or by t. Neglecting meson masses, one considers for reaction (| 104[) 
the light-cone vectors p\ and P2 as the vector meson momenta (2pi ■ p% = s). 

In this Sudakov light-cone basis, transverse Euclidian momenta are denoted with underlined letters. The 
virtual photon momentum q reads 

Q2 

q = Pi P2- (106) 

s 

The impact representation of the scattering amplitude for the reaction (|104p is [2(| [H| 

8+ioo 

M = J^I ^<^~^/ d ^rn\-U,-r_ + g)J ^(£f G .(fc,fc',r) (107) 



6- 



where G w is the 4-gluons Green function which obeys the BFKL equation [22j. G w reduces to 

within the Born approximation. We focus here on the 7* — > p impact factor $ of the subprocess 

5 (/ci,ei)+7*(<7) -> 5 (fc 2 ,e 2 )+p T (pi), (109) 

illustrated in FigJ^L, in the kinematical region where virtualities of the photon, Q 2 , and <— channel gluons kj_, 
are of the same order, Q 2 ~ kj_, and much larger than Aqcd- It is the integral of the S-matrix element 
£7 T 9^pt g w j-j- n reS p ec t to the Sudakov component of the t-channel k momentum along P2 . Closing the contour 
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of integration below shows that the impact factor can be equivalently be written (see FigOj) as the integral of 
the re-channel discontinuity of the S-matrix element S^ T 9 ~' PT 9 ■ 

+00 +00 

$ 7 ^ p (k, L-k) = e 7 " — \ —SZ 9 ^ p9 (k,r-k)=e~ f > i — / — Disc K S7 9 ^ P9 {k, r - k) , (110) 
2s J 2tt Is J In p 

-00 

where re = (q + ki) 2 denotes the Mandelstam variable s for the subprocess (|109|) . as illustrated in Fig[§h.. 
Note that the two reggeized gluons have so-called non-sense polarizations e\ — e\ = p2\/2/s . 
Considering the forward limit for simplicity, the gluon momenta reduce to 

/ Cl = LlL_ + - p 2 + fe X) fc 2 — — ^L^p 2 4. fc_ L) /j2 = fe| = fcj = — fc 2 • (HI) 

s s 

Finally, let us note that when writing (|llip we took an exact kinematics for the fraction of momentum along P2 ■ 
This kinematics naturally extends the usual Regge kinematics in the case where t-channel momentum transfer 
along P2 is allowed, which corresponds to the skewed kinematics which is typical of GPD studies. In usual 
computation of impact factors used in fc^-factorization, one usually makes the approximation that these two 
fractions are exactly opposite. Here we make such a choice in order to introduce skewness effects in a correct 
manner since these terms will contribute at the twist 3 order we are interested in. Note that within kx- 
factorization, the description of impact factor for produced hadron described within QCD collinear approach 
requires a modification of twist counting due to the off-shellness of the t— channel partons. Therefore, when 
here we say " up to twist 3" we only mean twist counting from the point of view of the collinear factorization of 
the produced p— meson, and not of the whole amplitude, e.g. 7* p — > pp or 7* 7* — > pp. 

In order to describe the collinear factorization of p-production inside the impact factor IjllOp . we note that 
the kinematics of the general approach discussed in section |n] is related to our present kinematics for the impact 
factor (|110p by setting p = pi, while a natural choice for n is obtained by setting n = uq = p%/{px ■ P2) (this 
latter choice for n, though natural, is somehow arbitrary as we discussed above in section Til Dp . 

We will now distinguish and make a comparative analysis of two different approaches: LCCF and CCF. We 
will show that these two results are actually fully equivalent to each other, when using the dictionary III El 



B. Calculation based on the Light-Cone Collinear Factorization approach 



1. 7I, —¥ pl transition 



Working within the LCCF, we first recall the calculation of the 7^ — > p L transition which receives contribution 
only from the diagrams with quark-antiquark correlators 6 . It is given by contributions from the p^ term of the 
correlators (fTOj) of the twist 2. Higher order corrections would start at twist 4, which is below our accuracy. The 



Hereafter, except for final results, we perform the computation for a meson which would be a one flavour quark— antiquark state. 
Its wave function is then restored at the very end. 
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computation of the corresponding impact factor is standard [23| . It involves the computation of the 6 diagrams 
of FigdJ The longitudinal polarization of the virtual photon reads, in the Sudakov basis, 

■t-gG-f+TrfV <112) 



while the momentum of the p reads 



and the longitudinal polarization of the p is 



o 

Tfl 

tf=rf + -rP$, ( 113 ) 



•s 



1 / „ ml 



(114) 

Consider for example the diagram (a) of Fig[U as illustrated in FigfTUl Computing the corresponding S'-matrix 
element, the corresponding contribution to the impact factor reads 

4s 2N c 2sJ J 2n [{ypi-q) z +ir)\[(k 2 + ypi) z +ti]\ 

o 



where the factor 1/4 is reminiscent from the Fierz identity, the factor 2/s comes from the normalization of the 

2 jV r 



non-sense polarizations. Finally, the color factor is due to the fact that when summing over the color of 



the t— channel gluons, the net color coefficient for 7*7* — * pp should be — l)/(4JVg) (due a Fierz factor 
l/N c when factorizing each of the two p DAs). Among the two propagators, only the second one, involving 
(&2 + VPi) 2 + if] — kx — k 2 x + irj has a pole in k, contributing when closing the contour below (therefore 
contributing to the discontinuity). The result is then easily obtained after extracting the corresponding residue. 
Diagram (c) provide the same contribution, since it can be obtained from (a) by the replacement x «-> x. 
Diagrams (b) and (d) vanishes for this twist 2 transition. Diagrams contributes only when closing the k contour 
above. Finally, the net result for the 7^ — > pl impact factor is, after taking into account the p° wave function 

S7£-/*( fc 2) = 5ab f dyVl {y) ^ ? . (116) 

V - ; V2Q 2N C J y ^ ny yyQ2 + k 2 V ' 



Note that diagrams (a) and (c) of FigH] are the only diagrams which contribute when computing the hard part 
by closing the k contour below. 



2. 7y — > pt transition 



We now concentrate on the 7^. — > pr transition, which impact factor will be one of the main results of this 
paper. The 2-parton contribution contains the terms arising from the diagram Figf?! where the quark-antiquark 
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k-2 + yi pi 



ViPi 



ki - y\ pi 




ki - yipi 



-yi pi 



fa 



Figure 11: The detailed structure of the diagram (bl). 



correlators have no transverse derivative, and from the diagrams Figj5l where the quark-antiquark correlators 
stand with a transverse derivatives. The computation of the diagrams of Fig0] for the 7^ — > pr transition 
goes along the same line as for the twist 2 7^ — » pl transition discussed above. The practical trick used for 
computing the contributions of FigjS] is the Ward identity 



d 



where 



1 



m — t/> — ie 



(117) 



where lines denotes fermionic propagators. This leads to an additional Feynman rule when inserting a derivative. 
The corresponding insertions are denoted with dashed lines in Figj5] 

Consider for example the diagram (bl) of FigQJ illustrated in FigQTj Computing the corresponding S'-matrix 
element for the vector part, the corresponding contribution to the impact factor reads 



12, 2 . S ab 1 



$ hi = - e <?T-(H)-9 fp m P 



4s 



2N r 2s 



dy 



dn Trtyj (#1 - yi>i) j* T (#1 - yii) foiifaih + vi>i)\ t 



2n 







[{fa - ypi) 2 + irj\ 2 [{k 2 + ypi) 2 + irj\ 



faiv), (118) 



This part of the impact factor receives (identical) contributions when closing the n integration contour either 
from above or from below, which reads 
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e q g 2 , S ab f ~e* T -e 7 {yyQ 2 + kf) + 2e* T -ke* T -k{l-2y) T 
— —f P m p — I dyy — — i2n9 <p 3 (y). 



{Q 2 yy + k 



(119) 



The computation of the axial part is similar (note the i/4 factor from Fierz) 



(f> A_ *2 f s 2, 5 ab 1 f, [d* Tr (h-vti)-y a (h-yti) fa fa 75 fa (fa + Ml a T( , 

^bl—~ e q 7 _ V l )9 Jp m P7TTr7r l a V I o ui = \2 1 ■ ~\2\ { 1 \2 1 ■ 1 6 e'pn ^3 \V) > 

As 2N c 2sJ J 2tt [{fa - ypi) z + ir]\ 2 [{k 2 + ypi) z + if)] 





and leads to 



(120) 



Zq9 
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fplTlp 



~2N~ r 



dyy- 



(Q 2 yy + k 2 ) 2 



(121) 



The contributions of 3-parton correlators are of two types, the first one being of " Abelian" type (without triple 
gluon vertex, see Figj6|) and the second involving non- Abelian coupling with one triple gluon vertex (see FigJTJ) 
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Figure 12: The detailed structure ol the diagram (aGl). 

or two (see FigJS]). Let us first consider the " Abelian" class. They involve two kind of Casimir invariants: 

l-T r (t c t a t b t c )=C F -^ ^Ca^-: (aGl), (cGl), (eGl), (fGl) (122) 
1 , / AT\ S ab 8 ab 

— Tr{tH a tH h )=\C F -^\—=C b — -. (bGl), (dGl), (aG2), (cG2), (bG2), (dG2), (eG2), (fG2), 

where the 1/N C comes from the Fierz coefficient when factorizing the quark-antiquark state in color space. 
Again, to illustrate the method, we consider the peculiar diagram (aGl) of FiglU illustrated in Fig[12] The 
vector contribution reads 

12,., 2 . s ab if,, fdK rr[^y (yi j*i - rf) ^ (wi j*i - *f) jfe (lb + ite^i) jfe j^i] 



2tt [{yipi - q) 2 + irj\[(y2Pi - q) 2 + irj\[{k 2 + 2/2 Pi) 2 + if)] 



x B( Vl ,y 2 ), (123) 
and equals 



Kgi = -^tt- fp m P WTF / <*J/i d V* ^T B (yi> ^) • (124) 



The corresponding axial contribution reads 

_ 12 2 5 ab 1 f dn Tr[^ (jnfr - j) la (y 2 ^ - j)f 2 (fe + y 2 j> x ) frfa] 

aG1 ~ 6 Us [l)9 Spmp 2N c 2 S J dyiam J 2ir{(y 1 p 1 -q) 2 +i V }{(y 2 p 1 -q) 2 + ir 1 }{(k 2 +y 2 p 1 ) 2 + lV ] 



x e a e , Tpn D( yi ,y 2 ), (125) 
and equals 

1 



Consider now the "non- Abelian" diagrams of Fig[7l involving a single triple gluon vertex. They involve two 
kind of color structure: 

2 N 1 fi ab 

" -{-i)Tr{t c t b t d )r d = -£ — — : (atGl), (dtGl), (etGl), (btG2), (ctG2), (ftG2) (127) 



N 2 - l v ' v ' J 2 C F 2N C 

2 N 1 f) ab 

~ -(- l)T r(t c t d t b )r d = -^-—— : (ctGl), (btGl), (ftGl), (atG2), (dtG2), (etG2) 



N 2 -V ' y ' J 2 C F 2N C 
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2/1 Pi - g 



fci + (yi -2/2) pi 




in vl 



Figure 13: The detailed structure of the "non-Abelian" (with one triple gluon vertex) diagram (atGl). 

where the 2/(7V c 2 -l) comes from the Fierz coefficient when factorizing the quark-antiquark gluon state in color 
space. Let us consider the diagram (atGl) of Fig|7l illustrated in Fig[l3] We denote as 



d v P{k) = g"P 



k v nP + kPn v 
k ■ n 



the numerator of the gluon propagator in axial gauge, and 

^Mi Ate pt3 (fa: ^2, fa) = (fa — fa)m 9fj. lf j. 2 + 



(128) 



(129) 



the momentum part of the 3-gluon vertex, where fcj are incoming, labeled in the counter-clockwise direction. 
The contribution of the diagram (atGl) then reads, for the vector DA, 



4> 



v _ 1 2 (-i)N c 2 



S ab If f dn 

~2N~Tsj dyi dy2 J 2^ Tr ^ J ( yi ^ 1 ~^ 7i/ (h + V2i>i)hi>i\ 



x 7? \2 i ■ TrTI — 71 \ — \2 1 ■ frTz — 7^ \2 , ■ 1 P2 e T 67%. 2/2) ■ (13U) 

lu/ipi - 3) + ^J[( fc i + (yi - vrivi) + mKfa + 2/2 pir + «?7j 

Note that for this diagram, as well as for all "non-Abelian" diagrams, one can easily check that only the g up 
part of (|128[) contributes. 

Closing the re contour above or below gives for the vector DA part of the diagram (atGl) the result 



atGl 



^ m " fp 2N c c; ' dyidy2 



(yi -2/2)2/2 



2/1 (2/1 k + 2/2 (2/2- 2/1) <2 2 ) 



er ' e-y B (Vi, 2/2) ■ 



(131) 



Similarly, the contribution of the diagram (atGl) reads, for the axial DA, 

1 



A _ * 2 H)jVc 2 , 

atGl— e 9 7"~77; 3 m pJp 



4 s 2C F 



<5 ab If f dn 

JN~2^ / dyidy2 2^ Tr ^ 7 ~ ^ 7l/ (^2+^2^1)^2^175] 



d vp (fa + (yi - 2/2)^1) V p \ a (-fa - (2/1 - y 2 )pi, fa, (yi - y 2 )pi) \ 



[{yiPi - q) 2 +ir)][(fa + (2/1 - V2)pi) 2 + irf[[{k 2 +2/2 Pi) 2 + irj\ 
and closing the re contour above or below gives 



P^e% pn D( yil y 2 ), (132) 



atGl 



(5 Qb N r 



dyi dy 2 



(2/1 -2/2)2/2 



2/i (?/i i + 2/2 (2/2- 2/1) <3 2 ) 



er ' e 7 Z%i, 2/2) ■ 



(133) 
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Figure 14: The detailed structure ol the diagram (gttGl). 



We consider now the "non-Abclian" diagrams of FigJSl involving two triple gluon vertices. They all involve the 
color structure 

9 IV f\ ab 

5 Tr\t c t d ]f cea f edb = — — . (134) 

N 2 -l 1 U J C F 2N C V ' 

For illustration, let us consider the diagram (gttGl) of FigJHJ illustrated in Fig[l3] It reads, for the vector DA, 

l 

1 2 (— i)N c S ab If f dn 

^ttai = -e«?7~ c C g2m P^2A^2^ / dyidy2 j 2^ Tr ^7(2/i^i-^)7^A]rf i/ ''(-'? + (l + ?/i-y2)pi) 

o 

Vp\ a {q - (1 + yi - 2/2 ki, -fc 2 + (2/1 - y2)pi) d a(3 {k 2 + (2/2 - yi)pi) 
[{yipi - q) 2 + iii\[{-q + (i + yi - 2/2)pi) 2 + ^][(fc 2 + (2/2 - y\)p\) 2 + £17] 
x v PTS {k 2 + (2/2 - yi)pi> ~ k 2, (2/1 - y2)v\)v\v T i B (yi, 2/2) ■ (135) 

It equals, when closing the k contour below on the single pole coming from the third propagator, 
s v e g g 2 . S ab N 1 f B(y u y 2 ) , 

^ugi =-— m rfpw- c c^Q2j d yi *V* ■ ^ ■ ( 136 ) 



The axial DA contribution from the diagram (gttGl) reads 

1 

$ ^tGi = ~ e g jj - 92mp ^ p JN~h J dyidy2 J ^ Tr ^'< ^ X ~ ^)7^^i75]d" p (-«+ (1 + 2/1 -2/2) Pi) 



x VpA«(g - (1 + 2/1 - 2/2)Pl, fcl, ~fc2 + (2/1 ~ 3/2 rf a/3 (fc2 + (2/2 ~ 2/l)Pl) 

[(2/1 Pi ~ q) 2 + -iv][(-Q + (1 + 2/i - V2)pi) 2 + irj\[(k 2 + (2/2 - 2/i)pi) 2 + w?] 
x V j 3Ta{k 2 + {y 2 -yi)pi,-k 2 ,(yi~y 2 )pi)p^ple\, Tpn B{y 1 ,y 2 ). (137) 

It equals, when closing the k contour below on the single pole coming from the third propagator, 

^9 2 . S ab N c 1 f D( yi , y 2 ) , 

W="— m pfp2N e C;^ Jdyidy 2 ^—e T -e,. (138) 
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All other diagrams of each class can be computed according to the previous examples. 

In order to present now the full result in a compact form, we decompose the result impact factor into spin- 
non-nip and spin- flip part. The non-flip part is proportional to 

T n .f. =-(e 7 -e* T ), (139) 

whereas the spin-flip part involves 

= (e 7 ■ fc ± )(4 ■ k±) (e 7 ■ e* T ) 
} ' k 2 2 [ ' 

We label the transverse polarizations as 

e(+> = =-±- [ei+le2 ] = -±- (o, 1, i, 0) , (141) 

e<-> = = -L [ ei - i e2 ] = -L (0, 1, -i, 0) . (142) 



They satisfy 



and 



In this basis, 



£ (A)* = £ (-A)_ (W3) 

e (+) e (+)* = e (+) c (-) =_i, and £ (-) e H* = e (-) e (+) = _i. (144) 



4 +) 4 +) * + & 4~ } * = 4 +) 4- } + 4 +) = -g± „ • (145) 

Decomposing the impact factor as the sum of spin-non-flip and spin-flip contributions 
and introducing the notations 

a = k 2 /Q 2 and C ab = - & ^"l" ^ T (147) 
one obtain the following results for the two bodies contribution 

^r, PT (k 2 ) = — — c F jd yi ( ( 2 ?/i - 1) vUsn) + 2 jft (i - yi) ^ (yi) + y£ (yi) 

_ 2 Cf 7 \ 



o 

2 a (a + 2 1/! (1 - yi)) ((2yi - 1) <p[ ( Vl ) + <p T A (yi)) 



yi (l - yi) (a + yi (1 - yi)) 2 



and 



l 
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o 

The three bodies contribution reads 



(148) 



<t>J^ {k 2 )= C^l CF j ^ 4a [{l _ 2yi)ip T [yi)+ip T M] (149) 

2 C F J (a+(l-yi) yi) 



™ / 3 C F J J \ a+ (1 - yi ) yi - y 2 + l)a+ y 1 (1 - y 2 ) yi a + yi (y 2 - yi ) / 

o o 

Vi$B( Vl ,y 2 ) f a (2 C f -iV c ) (22/1-1) | aN c (l- yi ) \ 

a + (1 - yi) yi V (yi - V2 + l) « + yi (l - y 2 ) y 2 « + yi (y 2 - yi) / 



+ (C 3 y S(yi,y 2 ) + C 3 A J D(yi,y 2 )) 



2 C F yi 1 



a + (1-2/1)2/1 i-yi 



2 



^Vc (1 - y 2 ) (yi - y 2 ) g : 
;i - yi) a + (1 - y 2 ) (y 2 - yi) 



+ C F + N c 
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and 



f - 3 y ~ ' 2 C F a + yi (1-yi) \a (1 



2C f 



x [<gD( yi ,y 2 ) (l + yi 



j/i - 2/2) + 2/1 (1 - J/2) ay 2 + yi (-yi + J/2) 

2/2) + C;T- B (2/i>2/2) (1 - yi -j/2)] ■ 



(151) 



The full result for the impact factor reads, after several simplications due to the use of the equation of motion 
and the symmetrical properties of 2 and 3-parton correlators, as 



er^ 2 ) 



-S- < "2 dy 




(a + 2yi (1 - 3/1)) a 
J/i (1 - J/i) (« + J/i(l-J/l)) 

A r>A.. .. M yi (1-J/l) " 



+2 / [C 3 V S( yi ,y 2 )-C 3 A i?(2/i,2/2 
A^c 1 



[(2i/i-l) + 

2 - JV C /C F 



C F y 2 a+ j/i (j/ 2 - 2/1) 
2/i 



"+ 2/i(l-2/i) 
- 2 / d Vl dy 2 [$B ( yi ,y 2 ) + C£d {y u y 2 
(2 - N C /C F ) Vl a 



a (j/i - 2/2 + 1) + 2/1 (1 - 2/2) 

2 + at c /c f 



2/1 



a + 2/1 (1 - 2/i) V" (2/1 - 2/2 + 1) + 2/i (1 - 2/2) 
JV C (j/x - j/ 2 ) (1 - y 2 ) 1 



and 



C p 



$}^ pT (fc 2 ) 



J/i 



4 / dj/! dj/2 




a + 2/1 (1 - 2/1 
(2 - iV c /C F ) 



a (1 - 2/i) + (2/2 ~ 2/i) (1 - 2/2) 



£ -5 [^( W ) - (2j/i - 1) ^(2/l)] 

(a + J/i(l- J/i)) 

- [Cg 4 ^ (j/x, j/ 2 ) (-j/x + 2/ 2 - 1) + C 3 y B (2/1,2/2) (2/1 + 2/2 - 1)] 
iVn 1 



(152) 



a (2/1 - 2/2 + 1) + 2/1 (1 - 2/2) C F y 2 a + yi (y 2 



2/iJ 



(153) 



The gauge invariance of the considered impact factor requires a special attention. The 7* — > p F impact factor 
is constructed in such a manner that it should vanish when k 2 = 0. This fact is a consequence of the gauge 
invariance of the impact factor. From our final formulas (|152|) and (|153|) . it is obvious to check that and 
$ n ./. indeed vanish when k 2 = since T n j. is a phase for flip transition, which is regular in the fc 2 — > limit. 
The vanishing of the "Abelian", i.e. proportional to C F part of (|152[) is particularly subtle since it appears as 
a consequence of the equations of motions ([551 ISO]) . Because of that some comments can be useful. Let us note 
that the sum of Eq. ([55]) multiplied by y\ and of Eq. ([50)1 multiplied by y\ takes the form 



22/12/1^3(2/1) + (2/1 - 2/iVi (2/1) + ^a(2/i) 
1 



(154) 



-2/1 / ^2^(2/1, 2/2) + C 3 A ^(2A, 2/2)] -J7l / dy 2 [-^B(y 2l y 1 ) + C^D(y 2 ,y 1 )] , 



from which, after integration over yi of the both sides of (|154p multiplied by we derive the equality 



1 

/ (2j/ij7iV3(j/i) + (j/i - 2/1)^(2/1) + <^a(2/i)) 
7 2/12/1 







1 1 
dyi 





dy 2 - [$B( yi ,y 2 )+t?D( yi ,y 2 )] • (155) 
2/1 



Now, by inspecting the expression (| 148[) in the limit «->0wc sec that only the first term in {...} survives and 
it has a form of the l.h.s of expression (|155|) . Similarly, by inspecting the expression (|150p in the limit a — > 
we see, that only the last line of this expression survives in the limit a — > and that the resulting expression 
coincides with the r.h.s of (|155p . Consequently, the non- vanishing terms cancel out due to the relation (|155() . 
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At the same time, the vanishing of non-Abelian (~ iV c ) part of (|150[) is the result of direct cancellation of the 
non-vanishing contributions of the diagrams (bGl), (dGl), (aG2), (cG2), (bG2), (dG2), (eG2), (fG2) of Figl 
(see Eq. (|122l) ) with the corresponding ones coming from diagrams of Fig|7] and Fig[5] containing triple-gluon 
vertices. Thus, the expression for the 7* — > px impact factor has finally a gauge-invariant form only provided 
the genuine twist 3 contributions have been taken into account, hidden in formula (|152|) when writing $n./. by 
the fact that we have used e.o.m., which explicitly relate 2 and 3 particles correlators. 

We end up this section with a comment about the problem of end-point singularities. Such singularities 
does not occur both in WW approximation and in full twist-3 order approximation. First, the flip contribution 
(|153[) obviously does not have any end-point singularity. The potential end-point singularity for the non-flip 
contribution (|152[) is spurious since ipj^x\), ip^xi) vanishes at X\ = 0, 1 (this is enough to justify the regularity 
of the result in the WW approximation), as well as B(x\,X2) and D(x\,X2). 



C. Calculation based on the Covariant Collinear Factorization 

We now calculate the impact factor using the CCF parametrization of Ref.0] for vector meson DAs. Let us 
outline basic ideas behind our calculation. We need to express the impact factor in terms of hard coefficient 
functions and soft parts parametrized by light-cone matrix elements. The standard technique here is an operator 
product expansion on the light cone, z 2 — > 0, which naturally gives the leading term in the power counting and 
leads to the described above factorized structure. Unfortunately we do not have an operator definition for an 
impact factor, and therefore, we have to rely in our actual calculation on the perturbation theory. The primary 
complication here is that z 2 — > limit of any single diagram is given in terms of light-cone matrix elements 
without any Wilson line insertion between the quark and gluon operators, like 

(V{pvM{z)i^mO) and (V(p v )\^(z)^A a (t z)^(0)\0) , 

we will call conventionally such objects as perturbative correlators. Actually we need to combine together 
contributions of quark-antiquark and quark-antiquark gluon diagrams in order to obtain a final gauge invariant 
result. 

One should stress that despite working in the axial gauge one can not neglect completely an effect coming 
from the Wilson lines since the two light-cone vectors z and n are not equal to each other and thus, generically, 
Wilson lines are not equal to unity. Nevertheless in the axial gauge the contribution of each additional parton 
costs one extra power of l/Q, therefore a calculation can be organized in a simple iterative manner expanding 
the Wilson line. At twist three level it is enough to consider the first two terms of such expansion 

1 

[z, 0] = 1 + ig J dt z a A a (zt) + 0{A 2 ) . (156) 


For instance, the quark-antiquark vector correlator can be written, for the general case z 2 ^ 0, as 

1 

(U(pv)|^(z) 7 ^(0)|0) = (V(pv)|$(*)7pMM0)|0)-i 5 J dt{V( Pv )mz) ltl z a A a {zt)^)\Q), (157) 



where we formally inserted the Wilson line in the r.h.s and performed its approximate subtraction according 
to (I156p . Then using relation (|4Tj) . we express the gluon field operator in the second term of (|157[) in terms of 
field strength, which gives us the (V(pv)\'>p(z)'y f j,G a p(t z)tp(0)\0) correlator. For the later we apply again the 
procedure of Wilson lines insertion (and its approximate subtraction) 

(V( P vM(zh^(0)\0) = (V(j,vM(z)^[z,0}^(0)\0) 

1 OO 

-igj J dtdae~ ea (V(p v )\ip(z)[z,zt + na}-f f _ l z a n' 3 G a p(zt + na)[zt + na,0}il}(0)\0} + ... , (158) 


where . . . stands for the correlators with more than one gluon field. 

Such correlator naturally appears (after Fierz decomposition) in the expression for the impact factor generated 
by the leading order diagrams of perturbation theory. Now we proceed to the extraction of leading 1/Q 
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asymptotic. It is achieved, due to the dimensional counting reasons, by the substitution of the off light-cone 
correlators by their light-cone limit where z t + n a oc z, z 2 — > 0. In this limit, using the CCF parametrization 
of section Hi B 21 for the light-cone correlators, one can deduce after some transformations that 

(p( Pp ) |^(z)7^(Q)|0) U^o = 

f p m p -ip^(e* ■z)fdye^P- z Hh(y)-h(y)) + e;fdye*y(P-^g i l\z) +... , 
o o 

where . . . stands for the contributions vanishing at twist 3 level, and 



(159) 



Hy) 



dat\ I da 2 



V(ai,a 2 ) 



(160) 



The physical polarization vector satisfies e • p p — (or e • p = since p p = p up to twist 3). On the other hand, 
the polarization vector of transversely polarized meson is chosen to be orthogonal to the light-cone vectors fixed 
by the external kinematics: e • tiq = 0. But one should take into account that the ex vector defined by (|3"Tj) has 
a non- vanishing scalar product with the vector no, 



e ■ z 



p ■ z 



(161) 



This relation was used to derive (| 1 59|) . 

Note that the zt+na cx z condition means actually that the vector z (which is an internal integration variable 
for the impact factor) is approaching during this limiting procedure the direction of the light-cone vector n, 
z oc n. One should mention, to avoid any misunderstanding, that it does not mean that we must put, say in 
Eqs. like (|159|) , the e • z scalar product equal to zero. What we actually do when performing the 1/Q power 
expansion is a Taylor expansion of scalar functions F(p ■ z, z 2 ), which depend generically on the two variables 
p ■ z and z 2 , with respect to the variable z 2 , whereas any scalar product of z with other vectors should remain 
intact. 

Performing a similar sequence of steps we obtain the following result for the axial-vector correlator at the 
twist 3 level 



lfp m f 



(/7(p p )|^(z) 7Al 75V(0)|0)| 2 2^ = 
mfh e* a p^z^ J dy e^- z ) (^ (y) - g { ° } (y)) + ^ aPl e* a p^n^p ■ z J dy jvb«)g<£> (y) 



(162) 



with 



y v 
9 ( ±\y)=4(£ J dat J da 2 



A(ax,a 2 ) 



(163) 



Comparing the obtained results (|159|) , (|162p for the perturbative correlators with initial parametrizations 
(|2gp we see that at twist 3-level the net effect of the Wilson line is just some renormalization of the h function in 
the case of vector correlator, whereas for the axial- vector we obtain in addition to the function renormalization 
a new Lorentz structure, the last term in (|162|1 . Nevertheless, we found that the last term in (|162|) produces at 
the end a zero contribution to impact factor. 

Let us now discuss gluonic diagrams which involves quark-antiquark gluon correlators, like 
(p(p p )\ip(z)jfj_A a (w)ip(0)\0) . Applying our procedure one can easily show that at twist 3 level 



( P (p p )$(z) 7ll gA a (w)ip(0)\0)\ 
{ P (p p )\^{z) lpl5 gA a {w)^{Q)\Q)l 



wocz,z 2 — *0 



= - m phpPn e *T a / D ® 



V(ai,a 2 ) 



,i(p-z)a 1 +i(p-w)a g 



(164) 



+0 = -im p f3P p ,e a f} lS n p p 



Da 



A(a 1 , a 2 ) i( v . z )a 1+ i( v . w ) c 
a a 



see Eq. (|42|43[) . Note that the first nontrivial effects induced by the Wilson line insertion start for such per- 
turbative correlators at the level of twist 4 only. Therefore taking into account such 3-partons contributions 
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is quite straightforward: one needs to calculate, projected in accordance with (|164j) . diagrams describing the 
production of collinear on-shcll quark-antiquark gluon state. 

One comment is here in order. Perturbative expansion generates, among others, diagrams where the gluon 
field is attached not to the internal part of the diagrams but to the "external" quark (or antiquark) lines. Such 
diagrams, in accordance to the logic of collinear factorization, should be factorized in terms of not 3-parton but 
2-parton correlators. Quasi-collinear gluon radiation appears at large distances and corresponding subprocess 
should be factorized not in the hard coefficient but is included in the soft part of the process, described by the 
2-parton quark-antiquark correlator. 

The internal variable z (and w) integration can be reduced to the Fourier integrals 

dtze'M = (2tt) 4 5 4 (0, J d*zz a e*^ = -i(2n)^S\l) 7 

where I stays here for some combination of the external and internal momenta 7 . The corresponding intermediate 
calculations do not contain principle difficulties both for the case of 2-partons and 3-partons contributions. Note 
that the contributions computed here have the same hard part as the one of the section IIII Bl except for 2- 
partons contributions with derivatives discussed in that section which have no counterpart here. In what follows, 
we simply give the final results. Then we will discuss in details an important issue related with the restoration 
of the gauge invariance for the final result. 
We use the notations 

k 2 N 2 



Q 2 ' 1 N 2 — I 
For the 3-parton contributions we obtain the result 



with 



and 



$«« 9 (a) = J Di 



4> 3 = - e 9 2 ™pfp i!i {$«ff»(o) + A$ 3 } , (165) 
V2Q 2 2N C 1 v ' ; v 1 



A$3 = _Z^ f {^ V (z 1 ,z 2 )(z 1 ^z 2 ) + CiA(z 1 ,z 2 )(z 1 +z 2 )} (166) 

l J Z\Z 2 Zg 



2a 



Z\Z 2 Z 2 



x |(C^(«i,« 2 ) +tfA(z 1 ,z i ))T n . f . 
+ (C3V( Zl , z 2 ) - CiA( Zl , z 2 )) 2 Zl T f . 




where Dz is defined according to (|39[) and where we have used symmetry properties of gluon distribution 
amplitudes under the exchange of quark momentum fractions: V{z\ 1 z 2 ) = —V(z 2 ,zi), A(zi,z 2 ) = A(z 2 ,zi). 
Let us discuss now the 2-parton contributions. We obtain 



$2 = _ eg m p f p S a , 2 . 



with 



<T»""(m)= / ,h { T n . f .9+(z) a _ (a + 2Z % + T f $-(z) I" , } (Ki-)) 
1 zz(a + zz) A (a + zz) 1 1 



The subsequent integration over the internal momenta with delta function derivative is done using integration by parts. 
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where 



$*(a) = (2z - 1) - ± g± (Z) . g± (Z) , (170) 



whereas for A$ 2 term we get 

$+(z) 



A&=T n . f .J dzUl\z) 



2zz 



(171) 



Note that $ 9 '(a) and $ ?l?ff (a) vanish in the limit a — > 0, whereas A$ 2 and A$ 3 do not depend on a. Now 
we need to demonstrate that A$ 2 and A$ 3 cancel each other. It guarantees the property of the impact factor, 
$(a = 0) = 0, which is directly related to the gauge invariance. 

One can now separate from A$ 2 the contribution A<I> 2 that is due to functions h(z) and g±(z), which 
originates in our method from the Wilson lines insertion procedure, 

A$ 2 = A$ 2 + A$ 2 , (172) 

where 

^ ' 

To calculate A$ 2 , it is convenient to present h(z) and g±\z) in the form 

l l 
h(u)=(V Jdt J Dz5{u- Zl -z g t) V ^ Z2) , g ( f\z) = 4(£ J dt J DzS(u - zi - z g t) A ^' ^ ■ 
o 9 o 9 

Using this relation one can easily found that 

2 ./ t Z g ^222 21^2 

The second term in (|172j) can be reduced to the form 

l 

- ^ / dz{ln(z). 9 Ti(z) + ln(z).g"( z )} . 
o 

where 



Then, we separate the WW and genuine twist 3 contributions to A$ 



6' 



A$ 2 = A$ 2W + A$^ 9e ", 



in accordance with 

gK(z)=g nww (z)+g Usm {z), g u {z) = g uww (z) + g U9en (z) . 
Using the explicit expressions for these functions in the WW limit 



1 2 

g nww [z)= f du h(uh g liww (z)= f du Mu y 
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one can easily found that the WW contribution to A$£, vanishes, 

Then, using results of Ref.0], which allow to express g^ gen (z) and g^ gen (z) in terms of 3-partons DAs, we 
found after some transformations that 

A $2gen = T n .f. f v V(z 1 ,Z 2 ) f gl ~ Z 2 _ ln £l£l J_\ + +A A ( Z l> Z 2) ( Zl + *2 _ ^^1^2 1 

2 J \ 3 z g \ ziz 2 z 2 z 2 z g ) 3 z g \ z x z 2 z x z 2 z g 

These results mean that 

A$ 2 = %^ / -J^{$V(z 1 ,z 2 )(z 1 -Z2) + tiA(z 1 ,z 2 )(zi + z 2 )} , (174) 

* J ZgZ\Z 2 

and thus that the constant terms of 2-parton (|174[) and 3- parton (|166|) contributions cancel each other 

A$ 2 + A$ 3 =0. 
Finally, the impact factor is given as a sum of two contributions 

3(a) = - 6g2 ™ p fp f Z {$™(a) + $«»( a )} (175) 

where $ 9 ^(a) and <& 9 ^ 9 (a) are given in eqs. (|169|) and (I1G7[) . 

D. Comparison of the two computations 

The above results for the 7* — * p impact factor were obtained based on the LCCF and CCF method, and look 
at first sight very different. As a testing ground of the validity of the dictionary elaborated in section Hi El it is 
interesting to show the exact equivalence between the two results. Since the detailed proof is rather involved 
and technical, this is done with whole details in Appendix [Bj 

IV. CONCLUSION 

In conclusion, we presented a general formalism which allows us to include in a systematic way higher twist 
effects. The general scheme was exemplified on the study of the r/10-meson production up to twist 3 accuracy. 
We did this using two different methods: the LCCF method in the momentum space and the CCF method in the 
coordinate space. The crucial point in the comparison of these two methods was the use of Lorentz invariance 
constraints formulated as the n-independence of the scattering amplitude within LCCF method, which leads to 
the necessity of taking into account the contribution of 3-parton correlators. We derived the dictionary between 
the DAs appearing in the two methods. Due to different calculational techniques used by the LCCF and CCF 
methods, we performed the calculation of the j*p impact factor in two different ways; we checked that they 
lead to the same result. This does not preclude the solution of the well known end-point singularity problem 
d, which may still require a separate treatment. The phenomenological application to our result to the 
HERA data [l(| will be discussed in a separate publication p5| . 
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Figure 15: 3-gluon vertex with one open index fi. 
APPENDICES 

Appendix A: n— INDEPENDENCE CONSTRAINT FOR THE N c STRUCTURE 

In this appendix we show that the n— independence constraint is automatically fulfilled for N c contributions. 
According to Eq.(|I52J, one should consider the "Abelian" diagrams (bGl), (dGl), (aG2), (cG2), (bG2), (dG2), 
(eG2), (fG2), and from Eq.([T2"7p and Eq.([I53]) all the non- Abelian diagrams. 

As a preliminary, let us consider a triple gluon vertex with a t— channel gluon of momentum k (which will 
stand for k\ or — k 2 in the all cases discussed below), saturated with a non-sense polarization p 2 , and a real 
gluon of momentum i g and polarization e g , leaving the index \i open, as illustrated in Fig |15l It reads 

V\ a (- k + Z 9 , k ' -£ 9 KP2 =-{2k-t g ) e g p% + (k + t g Up 2 e g ) -(21,- k)p 2 e£ (A.l) 

which reduces, after making the replacement e g — > £ g and using the fact that in our kinematics k has no 
component along pi , to 

Kai-k + ^g, k, -l g ) >t%I$ = -(2 h-Qp* - (t„ -p 2 ) [l g -kr = - [k 2 — (k — i g f] p» - {t g -p 2 ) {t g - k)» (A.2) 

where in the last line we used the fact that i g in on mass-shell. Hereafter, we will symbolically denote with the 
index p 2 the first term in the r.h.s of Eq. (|A.2[ ) and with W the second term (having in mind for this second 
term further use of Ward identities) . 

We now apply the same method based on the Ward identity which led to Eq. ([75)) . now for graphs (bGl) and 
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(bG2). This gives 1 




The P2 term obtained after applying the identity (|A.2j) to the diagrams (btGl) and (btG2) gives, respectively 




and 

yi - y-2 




(A.5) 



In the perturbative Regge limit we are considering here, the terms k 2 /(k — £ g ) 2 and k' 2 /(k' + £ g ) 2 can be 
neglected, since they are of the order of k 2 /s , and therefore Eqs. (|A.3|IA.4| lA~5j) exhibit an explicit cancellation. 



The signs in l.h.s of Eqs JA.31 IA.4I |A.4t are related to the signs in front of the corresponding N c coefficients, see Eqs. ll22lll27> . 
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Turning now to the case of graphs (dGl) and (dG2) one gets, using again the Ward identity, 




(A.6) 



The P2 term obtained after applying the identity (|A.2|) to the diagrams (dtGl) and (dtG2) gives, respectively 



yi - vi 



and 




Contributions (|A.6|) . (|A.7|) and (|A.8|) add to zero. 

The proofs for (aG2) goes along the same line and reads 




, (A.9) 
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The P2 term obtained after applying the identity (|A.2[) to the diagrams (atGl) and (atG2) gives, respectively 



and 




(A.10) 



(A.11) 



Contributions (|A.9[) , (|A. 10|) and (|A.11[) add to zero in the Regge limit. 
The proofs for (cG2) goes along the same line and reads 

yi -2/2 yi - V2 




(A.12) 



The P2 term obtained after applying the identity (|A.2|) to the diagrams (ctGl) and (ctG2) gives, respectively 



yi - y-2 



and 




Contributions (|A.12|) . (|A.13|) and (|A.14j) add to zero in the Regge limit. 
The proofs for (eG2) goes along the same line and reads 




(A.15) 
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The p2 term obtained after applying the identity (|A.2[) to the diagrams (etGl) and (etG2) gives, respectively 



and 




(A.16) 



(A.17) 



Contributions (|A. 15|) . (|A.16p and (|A.17p add to zero in the Regge limit. 
Finally, the proofs for (fG2) goes along the same line and reads 




(A.18) 



The P2 term obtained after applying the identity (|A.2j) to the diagrams (ftGl) and (ftG2) gives, respectively 




Contributions (|A.18|) . (|A.19|) and (|A.20j) add to zero in the Regge li mit. 

We now consider the "Ward" part (second term in the l.h.s of Eq. (|A.2[) ). first for the set (atGl), (dtGl) and 



41 

(ctG2). The contribution arising from (atGl) reads 




while the contribution arising from (dtGl) is 




which cancels the last term in r.h.s of Eq. (|A.2ip . The remaining term is thus the first term of the r.h.s of 
Eq. (|A.21j) . which equals the contribution arising from (etG2) since 
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Consider now the set (ctGl), (btGl) and (ftG2). The contribution arising from (ctGl) reads 



2/2 



-(V2 - yisvtJ, 



* r 


— - 

yi 






^ V2 — yi J 





- 2/2 



W 



£ g -P2 
{ki-Q 2 




yi °< yi 



yi - 2/2 



- 2/2 




yi - V2 
yi 




in 



- 1)2 



(A.24) 



while the contribution arising from (btGf) is 

-(y2-yi)Pix 



3/1 °\ yi 




-2/2 



Ig -P2 




V1 % Vl t g -p2 

(fcl - £g) 2 V^l - 2/2 _ „ , 



2/1 °L 2/1 



- y> 



W" 



yi - »/2 



(A.25) 

which cancels the first term in r.h.s of Eq. (|A.24[) . The remaining term is thus the last term of the r.h.s of 
Eq. (|A.24j) . which equals the contribution arising from (ftG2) since 



2/i - 2/2 



2/2 2/1 



2/2 , 




2/2 - 2/1 



-2/2 



2/2 2/1 



(fcl " ^ 




^7 = £ a- p2 
yi - 2/2 (fci - i g ) 2 




■w" 



- 2/2 



Consider now the set (etGl), (btG2) and (atG2). The contribution arising from (etGl) reads 
(V2 -Vi)p» 



2/2 ot yi 



-2/2 
(A.26) 




w 
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while the contribution arising from (btG2) is 



(2/2 -yi)P(i 



yi °L 2/1 

oJ 2/2-2/1 

~2/2 




(fc 2 +^) : 



r(^ + *a) A 




(A.28) 



which cancels the last term in r.h.s of Eq. (|A.27[) . The remaining term is thus the first term of the r.h.s of 
Eq. (|A.27|) . which equals the contribution arising from (atG2) since 



-(y2-yi)Pn 




°\ 1/2-2/1 




tg ■ P2 



- ffl -y 2 ifr+tgY 




2/1 



2/1 - 2/2 



- 2/2 



-2/1 



-2/1 



2/1 - 2/2 



Consider now the set (ftGl), (dtG2) and (ctG2). The contribution arising from (ftGl) reads 

2/1 

2/2 c 



(A.29) 



-(2/2 - yi)jv 




2 v-g 




tg -P2 



(fcl 



2/1 - 2/2 





(A.30) 



while the contribution arising from (dtG2) is 

^g ' P-2 




yi 



(kl ~ ig) 



2 (£ g klj/J 




Vl L-P2 



2/1 -2/2 ^s) 2 




2/1 H « 

2/1 - 2/2 



- 2/2 



M 2/2-2/1 



(A.31) 

which cancels the first term in r.h.s of Eq. (|A.30p . The remaining term is thus the last term of the r.h.s of 
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Figure 16: Two 3-gluon vertices connected, with one open index fi. 



Eq. (|A.30|) , which equals the contribution arising from (ctG2) since 
yi yi 



(2/2-2/1) P M 



V2, 




ft 

V2 - yi 



- y-2 



S/2 yi 



(k 2 + £ g y- 



"W" 




Collecting now the 4 remaining contributions (|A.23IA.26[ IA.29I IA.32|1 . one gets 



"Ward terms 11 = 2 



which vanishes since 



as expected from the fact that a quark-antiquark collinear pair cannot emit a gluon. 




(A.33) 



(A.34) 



We end up by considering the diagrams (gttGl), (httGl), (gttG2) and (httG2). As a preliminary, let us 
consider two triple gluon vertices connected by a gluon line, with t— channel gluons of momentum k and k' 
(with k = k\ (resp. k = —k 2 ) and k' = k 2 (resp. k' = — fci) for (httGl) and (gttGl) (resp. (httG2) and 
gttG2)) saturated with non-sense polarizations p 2 and with the gluon momentum £ g7 leaving the index [i open, 
as illustrated in FigfTHl It reads 

V^(-k+k'+£ g , k, -e g - k ')V\ a (£ g +k', -k', ~i g )e a g p x 2 p T 2 = -(2k-l g )e g p^+(k+e g )^{p 2 -e g )-(2£ g -k)-p 2 e^ 

(A.35) 

which reduces, after making the replacement e g 



£ g , to 



V^(-k + k' + £ g , k, -£ g - k') V\ a {l s + k', -k', -l g )t a g p x 2 p T 2 

= {t B ■ Pa) H (2k-k' + 2e g -(k-k')- (k' + £ g f) + (£ g ■ p 2 ){k' + £ g - k)»] (A.36) 

where in the last line we used the fact that £ g in on mass-shell. Hereafter and as for the case of a single gluon 
vertex, we will symbolically denote with the index p 2 the first term in the r.h.s of Eq. (|A.2|) and with W the 
second term (having in mind for this second term further use of Ward identities). 

Let us first consider the "Ward" terms. Applying the identity (|A.36[) to the diagram (gttGl) leads to 



(2/2 - 2/i) Pm 




- U2 



(VP2) 2 



(fci - k 2 



l {k 2 



(-fci + k 2 +£ g )^ 




,yi 



- U2 



"W" 



yi - yi 
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4G 



from which we deduce that the sum of the two contributions (|A.43[) and (|A.44|) equals zero, has was already 
shown in Eq. (|A.34|) . 

Applying again the identity (|A.36[) to the diagram (gttG2) leads to 




.VI 

- In 



i L V2-V1 



C2 P2fj, 




"p'i 



while the same identity applied to (httG2) gives 



.VI 



(v2 -yi)Pfi 




2/2 - yi 



-,V2 



C2 P2fj, 



"P2 



- y -2 




y-2 - yi 



V2 °L yi 



1/2 - yi 



- y-i 



(A.45) 



(A.46) 



from which we again deduce that the sum of the two contributions (|A.45|) and (|A.46|) equals zero. 

This achieves the proof that the n— independence constraint for the iV c terms is automatically fulfilled and 
does not lead to any new set of equations. 



Appendix B: COMPARISON BETWEEN LCCF AND CCF RESULTS 



Let us first consider the 2-parton spin- flip contribution (|149l) . Using the equation of motions IBT)]) we can 
express ipj and y T A in Eq. (|149|) in terms of the remaining distributions and we thus obtain from the LCCF 
approach 



/■2 



(B.l) 



dy\ 4 a 
(a+(l-t/i) yi) 



in 



((l-2y 1 ) 2 + l)<p 3 (y 1 )-(l-2y 1 ) l p A (y 1 ) + 2y 1 / dy 2 B - C? D)(y 2 , Vl ) 



Now our aim is to show that the Eq. (|B.1|1 coincides with the corresponding expression of the CCF approach, 
i.e. with the flip part of Eq. (|169p . which reads 



dz 2a 
(a + zz) 2 



z z 

'\gl + {£ J da! J da : 



(2z-l) 
A(a\,a2) 



dv($\\{z)-g\\{z))-<% / dai / da 



V(ai,a 2 ) 



(B.2) 



We begin with rewriting (|B.2[) using the results of Ref. Q in such a way that the form of resulting expression 
will be similar to one of the Eq. (jB.lj) . For that we reexpress in (|B.2I) the twist-2 distribution <?!>||(,z) in terms of 
other DAs. Let us note that </>|| (z) satisfies equations 



dz 



(B.3) 



or 



' (z) ^l^ )(2) -4^ " 



( a ) _ n l1gen 



(B.4) 
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which follow in an obvious way from the definitions of DAs and their solutions in the WW approximation. 
Eqs. (|B.3[ IB.4P lead in turn to the consistency condition 



(B.5) 



from which it follows that 



g i "\z) = 2g i l\z) + ~(z-z)^g i "\z)-2zg Ugen -2zg lUen +2 J dvg U9en (v)-2 J dvg" 9en (v). (B.6) 

o o 

This relies on the boundary conditions of the different DAs (see ) , which read 

g U {1) = ff iT (0 ) = g nWW {1) = g l1 WW {Q) = t g (a) (Q) = g (a ){1) = Q 



and 



l l l 

dug nww (u) = f dug n (u) = f dug^ ww {u) = I dug l ^(u) = 1. 





Note that they satisfy the symmetrical properties 

g^(u) g^ (u) . 



(B.7) 



(B.8) 



(B.9) 



The substitution in Eq. (|B.2|) of 4>\\(z) by the expression (|B.3|) and then replacing resulting g°^s by the formula 
(|B.6p leads to the equation 



*?*>)= T, 



C 3 A / dai / <ia 



(B.10) 

A(ai,a 2 ) 



o o 







in which we observe that the two terms in the first line of (jB.101) reproduce -using the vocabulary formulas (|1021 
11031) - the first two terms of expression in [...] of Eq. (|B. 1[) . 
Now consider the expression in the second line of (|B.10|I . i.e. 



2nd line 



dz 2a 
(a + zz) 2 



(B.ll) 



4, A A A 6 

2zV^«(z)-2z / <fe fl U«-»-tf(2*-l) [ da, [ da 2 + f ^ f doi2 ^ll^± 



Using results of Ref. the formula for 3>' 9 2 (a)| 2 „d u n e reads 



o o 



2nd line 2/ 



dz 4a 
(a + zz) 2 



(B.12) 



x C 3 V 



+C3 4 



' ± \z/± \ z z 

J * if («) + M(z) \- z J dvij ^ H{u) + M(«) J - i(2z - 1) /" rfaj / da 2 



\« 







V(ai,a2) 
(1 — ai - a 2 ) 2 



z 2 ( / * L(«) + A(z) J - zjdvlj^ L(u) + N(v) j + 1 / dai / da 2 



\u 







a 2 ) 

{I - ax- a 2 ) 2 
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and after interchanging of integration over the variables v and u it can be simplified and takes the form 



*f.(a)\ 2ndline =T f J 



(B.13) 



z 2 M{z)-z f dv(H(v) + M(v))--(2z-l) [ da x [ da 2 - V ( a ^ a2 \ 
J 2 J J (l-ai-a 2 y 





f 



z z N(z) - z dv (L(v) + N(v)) + - dan da 2 



A(a 1 ,a 2 ) 
' (1 - at - a 2 ) 2 



00 
The expressions M(v), H{v) 1 L(v) and N(v) are the integrals over the vector and axial- vector 3-partonic DAs 



V V 



M(v) = — f dai [ da 2 — 
dv J J 1 



V (ai,a 2 ) 
a.\ — a 2 







(B.14) 



V V 



da 2 ( d d , 

a 2 -i ycai- — V(ai,a 2 ) , 

— 0.1 — a 2 V " ,a 'i- " a i 







(B.15) 



V V 



N(v) = / da 



da 2 







1 — a.\ — a 2 \da 2 dai 



+ -j— ) A(ai,a 2 ) , 



(B.16) 



V V 



L{v) = ij " dai J t 



( a i3 a 23 — ) ^("l, "2) 

ai — a 2 V aai da 2 



(B.17) 





Using the vocabulary (|f 0Q[ If Of p and calculating the derivatives in above definitions we put them into the forms 

1 u 

M(u) = - / dvB(u,v)+ / dvB(v,u), (B.18) 



H(u) = I dv ( — u-^- ) B(u, v)- [ dv (— h u— ) B(v, u) , 



v — u du 



u — v du 



(B.19) 



N(u) = - / dvD{u,v)- / dv D(v, u) + 2 / dy 2 / dy 



u 



V2-yi 



(B.20) 



1 u /1 u \ 

f d f d if du f dv 

L(u) = —u I dv — — D(u, v) + U I dv —— D(v, u) — (u — u) \ I Z)(u, v) — I u) 

J du J du \ J v — u J u — v 



(B.21) 



from which we calculate integrals appearing in the formula (|B.13|) 



z J duM(u) = z J dyi \ dy 2 B(y ll y 2 ) 




(B.22) 
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z 

z J du H(u) = z 



J dyi J dy 2 B(y 2 , j/i) ( - 1 - + 1 j + : / du Lh ; . // ! + : / du. B{ u. : i 



(B.23) 



z / duN(u) = z / dyi / dy 2 D(y 2 ,yi) 1 



2(*-j/*) 
2/1 - 2/2 



(B.24) 



z j du L(u) = z 




z duD(z,u)-z duD(u,z) + / dyi / dy 2 D(y 2 , yi) -1 + 



2/2 ~ 2/2 

yi - 2/2 



2 

Finally, after substitution into Eq. (|B.13p of the expressions (|B.141 ... . IB.25|) and noting that 

z 1 



dai / da 



V(ai,a 2 ) 
(1 - at - a 2 y 



dyi / dy 2 



2/2 - 2/i 



(B.25) 



(B.26) 



and 



dai / da 2 



A{ai,a 2 ) 

(1 — ai - a 2 ) 
o o 

and the simplification resulting on the following properties 

1 2 l 



= dyi dy 2 



-P(2/i,2/2) 
2/2 - 2/i 



dzf sym (z) I dyi I dy 2 B ^ yi,V2 ^ = Q 



2/2 - 2/1 



(B.27) 



(B.28) 



and 

1 2 1 

rf2 fantisym(z) [ dyi [ dy 2 ^Ml^l. = 

J J 2/2-2/1 

2 

where f sym (resp. fantisym) is (anti)symmetric with respect to z <-> z , we obtain that 
$f»| 2 „ d / me =T / y (Q ^^ )2 [2^ I du (^B(u,z)~CiD(u,z)) 



(B.29) 



(B.30) 



which reproduces the last term in Eq. (jB.ip . 

We now proceed our comparison by considering the 3-parton flip contribution. We should thus compare the 
last line of (|167p . which is the covariant 3-parton flip result, with the corresponding LCCF result (|151[) . This 
is straightforward after using the dictionary (jlOOi II 1|) in the last line of (|167[) , in order to express it only in 
terms of B and D. Then, using the symmetrical properties Eq.(|2"7) (B(yi,y 2 ) and D(yi,y 2 ) are respectively 
antisymmetric and symmetric under the transformation yi <-> y 2 ), one should rewrite the prefactors of B in 
an antisymmetrical form (and correspondingly the prefactor of D in a symmetrical form). Doing this for both 
expressions leads to identical expressions. 

The same method applies to N c 3-parton non-flip contribution. Indeed, comparing the N c /Cp = 2cf term 
of Eq. (|150[) with the corresponding c/ contribution for the T n .f. tensor coefficient of (|167p after using the 
dictionary (jlOOl 1101} and restoring the symmetry of the integrands, one get identical expressions. 



We now focus on the non-trivial proof that the Cf contribution arising from the non-flip contributions are 
identical in LCCF and CCF approaches. Only this part leads to potential violations of gauge invariance, as we 
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saw in CCF approach of sectio n?!! CI where it was needed to combine A<& 2 and A<£> 3 . In LCCF approach, this 
also requires to combine 2 and 3-parton contributions. We should thus prove that the 2-parton contribution 
(|148p supplemented by the constant term (in color space) of Eq. (|150p is identical with the constant term (in 
color space) of the T n ,t, coefficient of (|167[) supplemented by the T„./. coefficient of (I169p . 

Consider first the 3-parton non-flip contribution with Cp color structure of the LCCF result l|150p . It reads 



^f f 7 TCF (k 2 ) = C ab Jd yi Jdy 2 
o o 
2a Vl (2t/!-l) 



(B.31) 



x |C 3 y 5 (2/1,2/2) 
+ ^D( yi ,y 2 ) 



1 



2 2/1 



2/2/1 + 1/2/1 



a +(1 -2/1)2/1 (2/1-2/2 + 1)^ + 2/1(1-2/2) " + (l-2/i)2/i 

-2ayi 1 | 2yi 2 1 

a + (1 - 2/i)2/i (2/1 - 2/2 + 1) a + 2/1 (1 - 2/2) a +(l-2/i)2/i Vi 2/i 



In this expression, we have written the last term — 1/j7i of B and D contributions as — 2/j7i + 1/j/i in order 
to prepare the identification with the CCF approach result. Removing for a moment the l/y~i contribution 
from Eq. (|B.3ip (i.e. last term for both B and D brackets), the obtained expression vanishes in the a — > 
limit (we therefore denote it with an index van.), and, after restoring the appropriate symmetry under 2/1 <-> y 2 
(antisymmetrized for B and symmetrized for D), it reads 



It ~P t c f { h 2 



n.f.3 van. 
1 



2/1 + 2/2 



1 1 

C ab a J dm J dy 2 |c 3 V S(2/i, 2/2) 


2/1-2/1 1 2/2-2/2 1 1 1 1 



(2/1 - 2/2) ((2/1 + J/2) a + 2/1 2/2) V2/2 



1 

2/1 



2/1 - 2/2 2/1 (" + 2/12/1) 2/1-2/22/2(0 + 2/22/2) 2/1 « + 2/12/1 2/2 a + 2/22/2 



-C3 4 D(2/i,2/ 2 ) 



1 



2/1 + 2/2 



a + 2/i2/i) V J/l(2/l ~2/2) 



(2/1 - 2/2)((2/i + 2/2)" + 2/1 2/2) 
1 1 

2/1 



1 

2/2 



1 

?7i 
1 



a + 1/22/2) V 2/2(2/1 - 2/2) 



1 

2/2 



(B.32) 



The corresponding expression obtained from the CCF approach reads, according to (|167p . 



qqg c F 

n.f. 1 



a) 



Dz 



Z\ Z2 Zg 



Z g Zi 



Z\z\ 



Z\Z 2 



az g + Z1Z2 zi(a + Z1Z1) a + z 2 z 2 



(B.33) 



which after using the dictionary (|100[ 1 1 1 [1 turns into 



4> 



r f CF (a) = -J dzi 


(zi + z 2 )z 1 



dz 2 



2a 



z\zi(\ - z\ — z 2 ) 

z\ z 2 



(C 3 v B(zi,l-z 2 ) + !-*&)) 



zi z\ 



a(zi + z 2 ) + z\ z 2 zi(a + zizi) a + z 2 z 2 



This reads, after restoring the symmetry properties of the integrands, 
1 1 



$™f*(a) = a / dyi / dy 



1 



2/1 - 2/2 



x C 3 / ^(2/i,2/2) 



2/1 + 2/2 



1 



(2/1 + 2/2) a + 2/1 2/2 V2/1 



+CiD{yi,y 2 ) 



2/1 + 2/2 



1 



(2/1 + yi)a + 2/1 2/2 V2/1 



1 

2/2 

1 

2/2 



2/2 



1 



2/2 



1 



2/1 



1 



2/1 



2/1 a + 2/12/1 
272 1 



(B.34) 



(B.35) 



(B.36) 



1 



2/2 a + 2/22/2 2/2 a + 2/22/2 2/1" + 2/i2/i 
2/2 1 , 2/i 1 , 2/i 1 



2/1 " + 2/1 2/1 2/2 a + 2/2 2/2 2/2 a + 2/22/2 2/i a + 2/i 2/1 
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which agrees with the result (|B.32|) after elementary algebra. 

We now prove that the 2-parton non-flip Cf LCCF result (|148|) agrees with the non-flip Cf CCF contribution 
®rif ' ai t er supplementing the LCCF result with the 3-parton l/yi term of Eq. (|B.31[) . 

Starting from the LCCF result (|148p and using the equations of motions (|59l [60)) we can express (pj and ip T A 
in Eq. (|149[) in terms of the remaining distributions, which leads to 



(jab 

~2 



dyi 
yi m 



2 a (a + 2 yiyi ) 



2 J/i2/i <P3(yi) ~yi dy 2 [(Y B(yi, y 2 ) + Ci D{y ll y 2 )} 



{a + j/i yi ) 

Vl J dy 2 [~C 3 V B(y 2 , Vl ) + C 3 A D{y 2 , Vl )] > - ViJ dy 2 B( yi ,y 2 ) + C? D{ Vl ,y 2 )] 



o 

Vl 



ill 



Vi / dy 2 B(y 2 , y x ) + C 3 A D{y 2 , j/i)] 



(B.37) 



Now let us consider the non-flip CCF result (|169i 1170]) , which reads 



K q f.(a)= T n.f / ^ 



+ J.9i Q) -C.3 4 / da, J da 2 



a{a + 2z z) 
z z(a + zz) 2 

z 

A{a 1: a 2 ) 



(2z - 1) 



dv 0||(z) - g { l ] (z)) - C 3 V / d "i / da 



V(ai,a 2 ) 



(B.38) 



The same series of transformations which led from (|B.38[) to (|B.10[) now gives 



KW= T n.f- / dz 



2a{a + 2zz) \ _ („) 



z — z 



C 3 / da\ / da 



z z(a + zz) 2 

V(oti, a 2 ) 1 A 



g ( l > (z)^zzg U9en {z) + z J dvg n9en (v) 

o 

A(ai,a 2 ) 



(B.39) 



C 3 / da\ I da 



Using the dictionary (|102[) one immediately gets agreement between the first term of (|B.37[) involving t/? 3 and 
the first term of (|B.39[) involving g\ . We now consider the remaining term of (|B.39[) . involving B and D through 

Using results of Ref. Q the formula for (a)\ rern . reads, after interchanging the integration over the 
variables v and u as we did in order to get (|B.13I) . 



K.fMrem.=T n . f . / dz 



2a{a + 2zz) 
z z(a + zz) 2 



(B.40) 



C 3 A 



1 



zzM(z) + z / dv(H(v) + M(v)) - -{2z - 1) / da x \ da 



V {0-1,012) 
{I -ax - a 2 ) 



-zzN{z) + z dv{L{v) + N{v)) / da x / da 



A{a 1 ,a 2 ) 
' (1 - ai - a 2 ) 2 



Finally, after substitution into Eq. (|B.40[) of the expressions (|B.14l ■■■ . IB.25[) and using the properties (|B.26[) 
and (|B.27|) . one obtains 

1 z 

$Z f .{a)\re m . = -T n j. J dz 2 -^±01 / dyi [g B(y u z) - D{ Vl , z)) . (B.41) 
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Coming back to the LCCF result, the remaining contribution after removing the ip^ term reads, after using the 
symmetrical properties (l27l) . 



-2 y dy 2 [-c 3 v yi) + C 3 A yi)] I • (B.42) 

J 

Now, the remaining term coming from the LCCF CV non-flip 3-parton contribution (|B.31|1 which we omitted 
for the moment is 

1 1 

<%7e^{k 2 ) =C ab J d Vl J dy 2 ^YB(y 1 ,y 2 )^- + CtD(y 1 ,y 2 )^ . (B.43) 







Using the symmetrical properties (|2T|) . this contribution is opposite to the contribution from the second line of 
(|B.42[) while the first line of (|B.42[) is identical to the result (|B.41|) ■ This ends the proof of the exact equivalence 
between LCCF and CCF results. 
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